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ABSTRACT 



Complementary binary sequences were invented by Golay in an in- 
vestigation of infra-red multi- slit spectrometry. This dissertation 
formalizes the basic results obtained by Golay and develops new concepts 
and techniques for examining the characteristics of these special binary 
codes. This work has developed new understanding of the structure and 
methods for the decomposition of complementary sequences. 

Complementary sequences have the property of an infinite correla- 
tion peak to ambiguity ratio when detected with a matched filter. These 
binary sequences should find much application as pseudo- random noise 
modulation signals for both radar and communications systems. 

A discussion of the need for such sequences is included in the intro- 
duction and is followed by a state of the art description. An operations 
group on the sequences is formulated and the proofs of several theorems 
concerning the operations group are given in a rigorous manner. One 
reason for developing the operations group is the application to elimina- 
tion of redundance in the computer search for new codes. 

Several invariant properties of complementary codes are proved 
through the use of the Hamming distance concept. Many more invariant 
properties of the sequences are demonstrated through the introduction of 
a Hamming vector. The concept of a Hamming vector is extremely use- 
ful as a complementary code decomposition tool. A large number of 
theorems are proved to enhance its use in this field. 

Several computer searches for complementary sequences are 
described and the actual computer programs for the GDC 1604 are in- 
cluded in the Appendix. 
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CHAPTER I 
INTRODUCTION 



Search radars of today are required to obtain smaller and smaller 
targets at ever increasing ranges* This has led to many different tech- 
niques of summing radio frequency energy to obtain these long ranges. 
Complementary sequences are well adapted for coding the RF energy 
pulses to give a summing of the returning RF energy in the echo. The 
present work is more concerned with the characteristics of complemen- 
tary codes such as their operation groups, Hamming distances and 
Hamming vectors, than it is with their application as modulating signals 
for RF carriers. However, the introduction will develop the radar 
problem as a vehicle to emphasize the use of these codes in future appli- 
cations in radar. Their application will probably be equally important 
in the communications field. The initial work in complementary sequen- 
ces was directed toward an infra-red application, but it is felt that the 
eventual application of these sequences will spread to many different 
fields. * 

The usual approach to the problem of increasing radar range is to 
look for ways of varying some of the parameters in the radar range 
equation which will lead to a greater maximum range. 

The radar range equation is: 



max 



GA 


SD 


t 


e 


P . 


1 6tt^ 


min 



1/4 



where = peak transmitted power in watts, 

P = minimum peak detectable signal in watts, 
min ^ ^ 

S = scattering cross section of target in units 
consistent with range, 

G = gain of transmitting antenna, 



A = effective area of receiving antenna in 
consistent units. 
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D = a composite loss factor for transmission lines, 
atmospheric losses, etc* 

The parameters under the control of the design engineer are A^,G, 

P . and A and G are fixed by various considerations such as 

min t e ^ 

antenna beamwidth, side lobe suppression, physical size of the antenna, 

etc. , and are not truly available as adjustable parameters beyond small 

variations. Therefore the two parameters for the designers to optimize 

are P . and P^. 
min t 

P . is basically limited by three factors: Johnson noise in the 

input circuits, shot effect and other noises in the first tube, and cosmic 

noise picked up by the antenna. P approaches KBT in a perfect 

min a 

receiver, where K is Boltzmann*s constant, B is the bandwidth and T 

a 

is the absolute space temperature, ^nrin currently being reduced 

toward the above value through the utilization of parametric amplifiers, 

masers and other low noise devices* Another method of reducing P 

° min 

is through pulse integration, since noise voltage averaged over a period 
of time has a value approaching zero while the signal values are additive. 
However, in most cases this integration of pulses can also be considered 
as an increase in transmitted energy. It will be considered from that 
viewpoint in this paper. There are two basic types of integration, co- 
herent and noncoherent. Coherent integration gives a gain of approxi- 
mately N, where N is the number of pulses integrated and noncoherent 

integration gives a gain of approximately a/N. Both of these values are 

3 

compared to a single pulse . The values are respectively the upper 
bound for coherent integration and the lower bound for noncoherent 
integration. 

There are basically three different integration time bases possible 
in a radar system. One is the antenna scan to antenna scan where the 
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operator or a computer, decides if the echoes from previous scans were 
signals or noise. The second method of integration is the interpulse 
method where the cathode ray tube, a delay line, or a computer sums up 
the hits within one antenna scan of the target. The third method is intra- 
pulse integration, or pulse compression, where different portions of 
one radar transmitted pulse are summed up. Complementary codes are 
used in one scheme of intrapulse integration. 

An examination of how the transmitted power or energy since 

A T, can be varied is now in order. From this relation, it is 
seen that there are two basic methods for increasing the transmitted 
energy: either increase peak pulse power, or increase the pulse width. 

A common method of increasing radiated energy has been to increase 
the peak power. This method has proven both expensive and wasteful as 
many radar transmitter tubes are peak power limited while only dissipa- 
ting a very small percentage of their average power capability, and 
radar modulators are taxed with a very similar problem. The alterna- 
tive did not look very promising at first, since to increase pulse duration 
normally reduces the high frequency content of a pulse, which in turn 
reduces range resolution capabilities. However, intrapulse modulation 
retains the high frequency content of a narrow pulse, while increasing 
the energy content to that of a pulse of long duration. Intrapulse modula- 
tion allows the summing, required in intrapulse integration, for 
detection in a radar receiver. 

There are two methods of intrapulse modulation, the digital method 
and the analog method. Although the emphasis in this paper is on 
digital modulation schemes - the analog system will now be briefly 
mentioned since it is a part of the overall picture. Figure 1. 1 contains 

4 

the graphs used in the explanation of the chirp radar scheme . In 
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Frequency modulation of a chirp radar* 
FIGURE 1. 1 
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Signal addition through the use of a delay line* 
FIGURE 1*2 
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this analog method of intrapulse modulation the frequency of the outgoing 



pulse is swept linearly upward with time as is shown in Figure I* 1-b. 
This gives a large high frequency content to the pulse while at the same 
time maintaining a long pulse duration. Upon reception the signal is 
passed through a reversed time delay, delaying the low frequencies more 
than the highs, as is shown in Figure 1. l-d« This allows all the fre- 
quency components of energy to arrive at the output simultaneously, 
giving a large and narrow output pulse. The type of filter used in the 
delaying process is known as a matched filter and is discussed in 
Appendix I. The procedure here is seen to be quite simple: first we 
generate energy components with various frequencies at different times, 
and then we delay all components so that they arrive at the output 
s imultaneou s ly. 

As a preliminary to the digital type intrapulse modulation consider 
a transmitted waveform such as in Figure 1. 2. This waveform con- 
sists of three pulses of unit length with spacing between them of one and 
three time units respectively. This waveform is then fed into a delay 
line in the receiver which has the delays shown in Figure 1.2-b. The 
output is either a one or a zero except at the exact match of signal to 
filter; the output is then three as is shown in Figure 1. 2-c. The three 
pulses have therefore been summed as though they were one, at the 
time between 6 and 7 units. These pulses still contain all the high 
frequency components of the short pulse of unit length but have effec- 
tively three times the energy transmitted in a single pulse at the time 
of exact match. There are two features to be noted about this output, 
first it has been delayed in time and second there are secondary pulses 
at various undesired times which might be classified as coding noise. 
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Fig, 1.3- Matched filter detection of a random code. 

Fig. 1.4 - Matched filter detection of an optimal noise code. 
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I will consider any undesired signal which might be mistaken for a true 
signal as noise. These two features will be discussed in some detail 
in the true digital intrapulse modulation schemes to be discussed later. 

The scheme just considered is somewhat inefficient since seven 
units of time are required to radiate three units of energy. If in the 
transmitted signal all the spaces with zero energy output in the pulse 
envelope had also had an output, seven units of energy would have been 
radiated in seven units of time. One might use some sort of frequency 
shift or phase modulation to accomplish this purpose, still utilizing of 
course the same output tube. Consider the same waveform as before 
with the units of energy coded as plus and minus by the modulator by 
means of phase shift modulation of either zero or 180^. The waveform 
is as shown in Figure 1, 3-a, Upon reception this waveform is now fed 
into a delay line in the receiver as shown in Figure 1.3-b. repre- 

sents a straight-through amplifier. Figure 1. 3-c shows the output of 
the delay line summer with respect to time. The output energy is now 
seven times that of a single pulse and has the frequency content, and 
therefore the range resolution, of one of the single pulses. A time delay 
in the output is again apparent, as are a number of noisy sub- peaks. 

This was a code picked at random and the results were good. There are 
certain classes of codes which are optimal, having even smaller noisy 

sub-peaks. These codes are of length 2^-1 and have been studied 
5 6 7 

extensively. * ^ They are known as L. codes or pseudo- random noise 

codes. An example of one of these pseudo- random noise codes is 
- + --+ + +, where the minuses and pluses indicate the phase of the RF 
signal. When this particular noise code is impressed upon its matched 
filter. Figure 1. 4-b, the output of this summed delay line is as shown 
in Figure 1, 4-c. It is to be noted that the maximum height of the 
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Matched filter detection of a complementary 
code of length 8. 

FIGURE 1.5 
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*'hash" (sub-peaks) is one, with a main peak of seven, whereas in the 
previous example picked at random the hash had a height of two. 

Elaspas has made a study and design of a radar utilizing this form of 

g 

modulation and matched filter detection. He infers that codes of 
reasonable length should be used. For example, a code of length 127 
would be a reasonable length to obtain the benefits of a pseudo- random 
matched filter radar. An optimum code of this length has a maximum 
subpeak of 13, which gives a good coding signal-to-noise ratio. 

There is a way, however, to completely eliminate the hash, al- 
though this required the transmission of two separate sequences. Both 
of these sequences could be transmitted through one output tube with a 

moderate band width. These codes which eliminate the hash are known 

9 

as complementary codes. An example of a complementary code pair 
being detected by a matched filter pair is shown in Figure 1. 5. The 
two codes being transmitted will be called A and B for convenience, 
with their respective matched filters also being labeled in a similar 
fashion. The receiver output, after the complementary codes are detec- 
ted by the matched filter and their outputs summed, is shown in Figure 
1.5-g. This demonstrates, as was mentioned before, that the hash level 
is zero. This lack of hash is indicative of an infinite coding signal-to- 
noise ratio. This particular feature is the prime advantage of complemen- 
tary sequence pairs over other forms of noise modulation. Figure 1. 5-g 

shows the same amount of delay as did all the other forms of matched 
filter detectors, and shows that the sum of all the energy sub-pulses is 
incorporated into the main pulse at exact match. 

Each of the systems described in this introduction has had the char- 
acteristic that the transmitted signal is spread over a wide frequency 
spectrum. This has in itself two advantages: One is that enemy detection 
of the signal becomes more difficult, and the other is that even after 
detection the jamming of a wide spectrum noise- like signal becomes 
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extremely difficult. Since the characteristics just described are de- 
sirable for a secure communication system as well as for a radar 
system, all of the advantages that accrue for the complementary se- 
quence modulation in radar would apply in a communications scheme. 

The initial major objectives of this investigation were two in 
number. 

1. The formulation of a set of operations on complementary 
sequences which form a group, along with a theoretical 
study of the invariant properties of such a group. 

2. The search for new complementary codes* An exhaustive 

computer search for a new code of length 26 was made 
as well as a partial search for codes of length 34. 

During the investigation two more major objectives were added 
to the work. 

These were : 

3. The application of Hamming distance to both the comple- 
mentary sequences and their group formulation^ 

4. The application of Hamming vectors to complementary 

sequence pairs and also to their group formulation^ 

The original objectives were oriented toward the generation of 
new codes, while the last two objectives were oriented toward the 
decomposition of codes* All of these objectives were attained and 
consititute a large portion of the rest of this work* The next chapter 
is a study of complementary sequences designed to bring the reader 
to the state of the art, and contains some proofs which have not ap- 
peared e Isewhe Chapte r 3 is the fo rmulation of the 
operations group on complementary sequences and, along with several 
theorems and proofs it also contains the identification of this operations 
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group with a known group« Chapter 3 as written requires only a min- 
imum of group theory knowledge* Chapters 4 and 5 have to do with 
the Hamming distance and Hamming vectors of complementary sequen- 
ces both in standard form and in group form. Chapter 6 is a short 
chapter and is concerned with the proof of two theorems used in the 
searching for new codes* Chapters 7 and 8 describe the computer 
search for new kernels of length 26 and 34 respectively, the actual 
computer programs used being listed in the appendicies. The last 
chapter, Number 9> contains the summary, conclusions, and suggestions 
for future research in this interesting field* 
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CHAPTER II 

COMPLEMENTARY SEQUENCES 
This chapter is presented to provide a background in complemen- 
tary sequences to enable the reader to better understand the material 
in later chapters. Most of the materia] in this chapter is an adaptation 
of Golay*s **Complementar y Series** with expansions and deletions to 

9 

fit the needs of this paper. The formalization of some proofs and the 
addition of new proofs given here have not appeared in print elsewhere; 
however, it is obvious from the tone of Golay*s work that he was aware 
of these proofs. This modification of Golay*s work is offered for com- 
pleteness of the dissertation and not as new work. 

In the introduction a pair of codes of length n were discussed which 
had the characteristic that when detected by matched filters, the sum of 
the two filter outputs was everywhere zero, except where the signal 
patterns had zero time delay with reference to the filter patterns. At 
the zero pattern delay time the output was Zn, leading to an infinite 
coding signal-to-noise ratio. A study of the necessary and sufficient 
conditions for this to be true will now be considered. 

Assume a complementary code pair (Aj,B) and suppose that one of 
the pair, i. e. A, is longer than the other. If this be true then there will 
be an output from the A filter with no cancelling output from the B 
filter. This is not allowed; therefore the A and B codes must be of 
equal length. Therefore assume a code A, n bits in length, with 
elements a^ through a^, where each a^ is either +1 or -1. Also, 
assume a second code B, n bits in length, with elements b^ similarly 
defined. Detect each of these codes with a matched filter (the time 
inverse of the code) and then sum the outputs of the A filter and the 
B filter. At time 1 (see Figure 2. 1) the signals a^^ and b^ have 
entered the filter segments a and b respectively. Since at this 



13 



a . . . 

n 



V-- 



Time 



Time 



Time 



. a^a^aj 



time 



•Vzh 




The a, and b. of the filter are 
1 1 

either ± I indicating straight thru 
amplifiers or inverterso 



1 a^a3a^a^ 



a a 



signal 



Ct 1 CL ^ 

n n- 1 n- 2 



a, filter 



bsb^bi 



b b 



signal 



n n- 1 n- 2. 



filter 



^ 4 ^ 3 ^ 2^1 



a a 1 ^ -> 

n n-1 n-Z« 



'>4'’3‘'2'=l 



b b , b ^ 
n n-1 n- 2 



r a 



n 



^n+l-r ^n+l-(r4l) 



Matched filter detection of a complementary 
pair of length n^ 

FIGURE 2.1 
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time the delay between signal pattern and filter pattern is not zero, the 

summed output voltage must be zero. This gives the condition + 

b^b^ = 0. At time 2 (Figure 2. 1) the signals are one segment farther 

into the filter and the output must again be zero, giving 

a,, a + b, b , + b 9 b = 0. 

2 n 1 n- 1 ^ n 

At time 3: a, a 9 +a.,a ,+a^a +b,b ^+b^b i+b^b = 0 

1 n-2 2 n-1 3 n 1 n-2 2 n- 1 3 n 

At time r ( r ^ n) : a ^ a , , +a^ 3* , i i \ . 

1 n+l-r 2 n+l-(r-l) . +a a +bib ,, +» o . 

r n 1 n+l-r 

+ b b =0 



r n 

At time n: a, a, a, + .«..o« + a a +bTbi +o*»«««+b b =2n 

nni l nn 

for at time n the pattern of the signal must exactly match the pattern 
of the filter. The explicit statement of the conditions given above for all 
times is 



F. = 
J 



i=n-j i=n-j 

E a.a^ ^ E b,b..^, 
1 = 1 ■’ 1=1 



i i+j 



= “ ( 2 . 1 ) 

= 2n j = 0 



An examination of the equation at time 1 (j= n - 1) shows the following 
eight possible solutions: 



ai 


a 

n 






1 


1 


1 


- 1 


1 


1 


- 1 


1 


- 1 


-1 


1 


- 1 


- 1 


-1 


- 1 


1 


- 1 


1 


1 


1 


1 


-1 


1 


1 


- 1 


1 


- 1 


- 1 


1 


-1 


-1 


- 1 



These possibilities show that if the A pair are alike » the B pair must 
be unlike, and if the A pair are unlike, the B pair must be alike. The 
equation at time 2(j=n-2) indicates that if the number of likes (like pairs, 
i.e. a^^and a^ p) of the A code and the number of unlikes (unlike 
pairs) of the B code are equal, and if the number of unlikes of A 



equals the number of likes of B, the equation is satisfied. Extending 
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this reasoning to time r shows again that if the likes of A equal the 
utilikes of B and if the unlikes of A equal the Jikes of B the equation 
is again satisfied^ 

It will now be proved that the last half of the previous statement is 
not necessary, since the first half of the statement contains the necessary 
and sufficient requirement. 

Let = number of unlike pairs in A for the spacing specified. 

(j in equation 2„ 1) 

Let L = number of like pairs in A for the spacing specified^ 

(j in equation 2. 1) 

Let = number of unlike pairs in B for the spacing specified* 

(j in equation 2„ 1) 

Let Lj^ = number of like pairs in B for the spacing specified. 

(j in equation 2'* 1) 

Assume that L^ = U^, we want to prove that = 1^^ 

Since the sum of the number of like pairs at each spacing with the num- 
ber of unlike pairs at the same spacing must equal the total number of 
pos sibilities , 

+ L^ = n-(n-r) for r ^ n where n is the length of the code 

and (n-r) is the number of units of time delay 
from exact match* 



For the same reason 

= n-(n-r). 

But since U, = L , 
b a 

U + U, = n-(n-r)o 
a b ^ ^ 

Therefore U = L, . 

a b 

The definition of a complementary pair of sequences will now be 

given. A pair of binary sequences of equal length with the number of 

like pairs of one sequence equal to the number of unlike pairs of the 

other sequence for each possible spacing is said to be a complementary 
code. 
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Two schemes of representation will be used in this paper* These 
schemes will be used interchangeably at the convenience of the author. 
The first representation is the one already presented where each a^ 
and b^ is either +1 or -1, The operation used in this case is 
ordinary multiplication. The second representation will use the opera- 
tion of modulo two addition denoted by where the elements a^ and 
b^ are either 0 or I, That these two groups (1 , - ) and (0,1; ® ) 

are insomorphic is shown in the following tables. 



• 


1 


- 1 


® 


0 


1 


1 


1 


- 1 


0 


0 


1 


-1 


-1 


1 


1 


1 


0 



where I -* 0 

-1 1 



The first form will normally be used when correlating of the codes 
with their matched filter is under discussion while the second form will 
normally be used when discussing the intrinsic properties of the indi- 
vidual codes. 

The necessary and sufficient condition for a pair of sequences to be 

complementary, is in terms of the modulo two representation 
i=n-j i=n-j 

F. = S + for all j, l<j<n-l (2.2) 

or an equivalent formulation by change of summation index 



f . 

J 



i=J 

S 



s 



. , (a.® a , (b.® b ... ® 1) 

1=1 ' 1 n-j+i' 1=1 ^ 1 n-j+i ' 



(2.3) 



The above S is standard summation. If modulo two summation is 
desired the symbol E will be utilized. (The equivalent conditions are 
given since they are both used in formal proofs in a later portion of 

the paper. Due to the possibility of confusion it seems more appropriate 
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to introduce the second form now rather than later when needed* ) 



U = L<i has been used to represent these same equations, and will be 
a D 

used elsewhere in the paper when it is more convenient to use it. 

An interesting property of any system of codes is the possible in» 
variances it may have under various types of transformations* It 
therefore seems appropriate to consider the possible transformations 
on a complementary sequence pair which leave the padr complementary 
and of the same length after the transformation* A necessary tool for 
this study is the property that a^fa^^^ = a^ f ^i+j symbol 

a. means the complement of a^, 0 = I> 1 = 0* 

Theorem 2, I 

The modulo two sum of a pair of binary numbers is equal to 
the modulo two sum of their complements* 

The truth of the theorem is obvious from the following exhaustive 
table of four possibilities: 



a. a. . . 

1 i+j 

I ® 1 = 0 

I ® 0 = I 

0 ® 1 = 1 

0 ® 0 = 0 



Therefore a. ®a. , . = a. 

1 i+j 1 



a. 



a. w a. , . 

1 i+j 

I«T = 0«0 = 0 

I ® 0 = 0 ® I = I 

0 ® I = 1 ® 0 = I 

0 ® 0 = 1 « I » 0 



i+j 



In each of the following proofs it will be considered that A and B 
are a complementary sequence pair of length n* Complementing each 
element in a code is called complementing the code* 

Theorem 2* 2 

Complementing the A code, or the B code, or both the 
A and B codes, results in a pair of complementary codes* (The 
proof is given just for the A code but is identical for the B code and 
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from both of these for the A and B codes, ) 



By hypothesis 
i=n-j 

F. = S (a.® a.^.) = 
J 1 i+J 

i= 1 



i=n-j 

S (b «b,^. 

1 i+J 

i=l 



1 ). 



Complementing the A code transforms 



a. - a. , 
1 1 



but by Theorem 2, 1 



a. ® a. . . = a. ® a. . . . 

1 i+j 1 i+j 

Upon substitution 
i=n- j 



i=n-j 



F. = S (a. «a.^.) = ^ (b. 

^ ‘ i=l ‘ 

Therefore complementing the A code has no effect on complementarity, 

-i 

and similarly for the B code or for both codes* 

>!C >1C 

Interchanging the first and last bits of a code, the second and next 
to last bits, the third and third from last bits, and so on is called time 
inversing or time reversing a code* 

Theorem 2* 3 

Time inversing the A code, or the B code, or both the A 

and B codes, results in a complementary pair* (This proof is again 

just for the A code, but also applies to the B code and therefore to 

both the A and B codes* ) 

By hypothesis 

i=n-j i^n-j 

f. = S (a. « a.^.) = ^ (b « b.^. ® 1), 

J i=i 1 i+J i=l ^ 



The time inverse of A causes each 

a . a . I . c 
1 n+l-i 

Applying this time reversal just to the A code shows that 
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i=n-j i=n-j 

F. = S (a. ®a.,.) -• I! (a 
1 i+j 



n+l-i ^ ^n+l-(i+j) 



J i=l " i=l 

Expanding both expressions for a few terms and remembering that 
modulo 2 addition is commutative as is regular addition we see that 
(on the left) (on the right) 



^n-1 


= aj 


® a 


^n-1 


= a ® 


a 


n 


n 




^n.2 = 




^n-1 +^2® 


F , 


= a ® 


a^ 


n- 2 


n 


t 



a , f a, . 
n- 1 1 

These expressions indicate that by a change in summation index the two 
forms of Fj are equal, as is the case* Therefore time reversing the 
A code has no effect on the complementary property, and similarly for 
the B code and for both the A and B codes^ 









The altering of a code is a transformation where every other bit of 
both codes is complemented. In a different portion of the paper a dis- 
tinction is made between altering odd bits and altering even bits. In the 
following proof, however, just the even bits are altered, although a 
slightly modified proof would hold for altering the odd bits. 

Theorem 2. 4 

The result of altering a complementary pair of sequences is 
again a complementary sequence pair. 

Since the codes are a complementary pair, 
i=n-j i»n-j 



F. = S (a. ® a.^.) 

j i=i 1 i+j' 



i; (b. ® b. . . ® 1). 
i=l " 



Divide the F. into two groups, F, for j= an even number and F 
j K r 

for j= an odd number. The portion of the above equations for j= an 



even number is 
i=n- k 



i=n-k 



1=1 1=1 
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For the theorem to be true, when j is an even number the above set 
of equations must be satisfied for all k when the pair is altered* 

(The proof is given only in terms of A, but duplicate steps must also 
be applied to B, ) 

When i=an even number, then i+k=an even number* In the alter- 
ing process since all even bits are complemented a.® a.,, -*a.®a.,, , 
but by Theorem 2* 1 these are equal* When i=an odd number, then 
i+k = an odd number and these are not affected by the altering process, 
since only even bits are complemented* Therefore F is invariant 

tv 

under altering, or symbolically Fj^ = ^kalt" 

Examination of F^ where r is an odd value of j, shows that of 
the pair a., one of the bits is complemented and the other is un- 

changed. 



a. ® a. , (a.® a. . or a.® a. . ) = a. ® a, . 

1 i+r 1 i + r 1 i+r i i+r 

Therefore F _ = (n-l)-F for both the A code and the B code. 

This signifies the likes of the A code for odd spacings are changed into 

unlikes, and the unlikes of the A code for odd spacings are changed into 

likes; however, for each chamge in A there is an opposite change in 

B so the total likes in A still equals the total unlikes in B for each 

spacing* It is to be noted that altering both the even and the odd bits is 

the same as complementing both codes* 












The operation of interchanging the A and B codes, although 
trivial, is also a transformation* There are a total of sixty four pos- 
sible transformations* These can all be generated by combining the 
above listed transformations* They will be discussed much more thor- 
oughly in Chapter 3* 
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Whenever a new facet of a pair of sequences such as the comple- 
mentary property is discovered, it is natural to wonder just how general 
this property might be. One question of interest would be the possible 
length limitations of the code pairs that might be forced by the comple- 
mentary property. The following proof due to Golay shows that 

2 2 

complementary sequences must have a length n of the form n=x +y , 

where n is an even integer and x and y are integers. The first 

theorem to be proved demonstrates that n must be an even integer 

2 2 

and the second theorem shows the x +y form. 

Theorem 2 , 5 

A necessary condition for a sequence pair to be complementary 
is that their length be an even number. 

By hypothesis 



i=n-j 

F. = E a.§a. , . 
j i=i 1 



i=n-j 

E b.«b.^. » 1 
i=l ‘ 



Let G. = F. modulo 2, so that 

^ i=n-j i=n-j 

G. = S a.® a.. . = S b. ®b. , . ® 1 . 

J 1 i+J 1 i+J 

An expansion of G for a few different values of j will suggest the 

proof. Start with j=n-l. 

G , = a,®a =b,®b ®1, 
n- 1 In In 

G ^ = a,®a ,®a,®a = b,®b ,®l®b,®b ®1, 

n-2 ln~12n ln-1 2n 

G , =a,®a ^®a , ®a i®a,®a =b,®b ,®1 ®b ,®l®b-®b ® 1 . 

n-3 1 n-2 n-1 n-1 3 n 1 n-2 n- 1 3 n 

Change the order of the sums in G^ 2 and G^ ^ to achieve the following 

2 



G -,=ai ®a^®a ,< 

n-2 1 2 n-1 



i =b, 
n 1 




g 1. 

1 



G -=a,®a,®a,®a ,®a , ®a =b,®b.®b,®b^ .®b^ ,®b fi 

n-3 1 2 3 n-2 n-1 n 1 2 3 n-2 n- 1 n 

In general for j=n-r 



G — ’a.®a-j®o,«®a ®a . .-i % ®a ® 1 ®b^®. « 

n-r 12 r n-(r + l) n-(r+2) n- 1 n 1 2 
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„..®b ®b ,..,®b - ..®b ,®b ® E 1 , 

r n-(4 + l n-(r+2) n- 1 n \ ’ 

and for j=n-(r + l) 

^n-(r + l)^^l®^2®- •” ®^r®^r + l®^n-r®^n-(r+l)®* • ' ®^n- 1 ®^n"^l ®^2^ 



r+1 



..®b ®b .,®b ®b , ^ ,®b - ...® 

r r + 1 n-r n-(r+z) n-(r+Z) 

Forming the sum (r+1) leave only five terms in the 

equation since G and G / diffe r only in a , a b .t, 

^ n-r n-(r + l) ^ n-r* r + 1" r+1" 

and a, bcause 
r r+1 

Si® S 1 = 1 . 

1 1 

G ®G / L \ - 3i = b , 1 ®b ® 1 , 

n-r n-(r+a) r+1 n-r r+1 n-r 



Adding ^)^^j®b^ ^ modulo Z to both sides of the equation give s 

a ,,®a ®b ,,®b =d. (2.4) 

r+1 n-r r+1 n-r \ • / 

Let n=Zs-l and r = s-l; then to satisfy the equation just derived 



a ®a -5 1 , , ®b ®b -5 , , , = a ®a ®b ®b = 1 , which 

s Zs-l-s+1 s Zs-l-s+1 s s s s 

is obviously false. 

Therefore n cannot have the form Zs- 1 and must be an even number. 

5jc >:< 



Equation Z^ 4 which was a step in the previous proof, is a very 

important necessary condition for a pair of codes to be complementary 

and will be referred to in the future as the parity check. The next 

Z Z 

theorem which shows the n=x +y form, starts with the basic assump- 
tion that A and B are a complementary pair of length n. 

Theorem Z. 6 

A necessary condition for a pair of codes to be complementary 
is that their length be the sum of the squares of two integers. 

By hypothesis 



i=n-j i-n-j 

F. = E (a. ®a.^.®l) = E 

j i=i 1 ^+J i=i 

Assume that the A code satisfies F. 

J 



and that the A sequence contains 
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ly 




I 




p ones and n-p zeros* The A code has a weight of p* Similarly 

assume that the B code satisfies and that its sequence contains 

q ones and n-q zeros* The B code has a weight of q. For the total 

likes in the A code, each bit is matched against each other bit exactly 

once. For the total unlikes in the B code each bit is again matched 
against each other bit once. There are two possibilities of likes in the 
A code, 1 matched against 1, and 0 matched against 0, The total 
number of like pairs for the A code is therefore the number of com- 
binations of p ones taken 2 at a time plus the number of combinations 
of n-p zeros taken 2 at a time, for the unlike pairs in the B code, 
q ones are matched against n-q zeros. 

This gives the equation 

PlPlil + <HlPH?:P-li = q(n-q) 

expanding and simplifying 

n=n^- 2qn+2q^+2p^- 2np, 

adding and subtracting 2pq gives 
2 2 2 

n=n -2qn+2q +2p - 2np+2qp- 2qp, 
collecting terms gives 

n= n^+q^+p^ +2qp-2np-2nq + p^ -^pq +q^ 
combining terms gives 
n= (n-p-q)^ + (p-q)^ 

Therefore complementary sequence lengths are permissible only in 
lengths which are even numbers and formed by the sum of two squared 
integers, 

5 *^ ❖ ^ 

A list of all possible code lengths up to 200 with the number of 
ones allowed in the A and B codes is given in Table 2* 1, 
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*Indicates possible kernel 



Code 


Uno rdered 


Code 


Unordered 


Lengths 


Weights of(A,B) 


Lengths 


Weights of (A, B) 


*2 


(2, 1)0.0) 


104 


(58, 56)(58,48)(48,46)(56,46) 


4 


(3,3)(1, 1)(3, 1) 


*106 


(60, 55)(51 ,46)(60,51)(55,46) 


8 


(6,4)(4,2) 


1 16 


(65, 6l)(55,51)(65,55)(6l,51) 


*10 


(7,6)(4,3)(7,4)(6,3) 


*122 


(67,66)(56,55)(67,56)(66,55) 


16 


(10, 10)(6, 6)(10, 6) 


128 


(72, 64)(64,56) 


*18 


(12.9)(9,6) 


*130 


(72,69)(72,61)(61,58)(64,57) 


20 


(13, 11)(9, 7)(11, 7) 




(73,66)(73,64)(64,57)(66, 57) 


*26 


(16, 15)(11, 10) 


1 36 


(76, 70)(66, 60)(76,66)(70, 60) 




(15, 10)(16, 11) 


144 


(78, 78)(66,66)(78,66) 


32 


(20, 16)(16, 12) 


*146 


(81, 76)(70,65)(81,70)(76,65) 


*34 


(21, 18)(16, 13) 


148 


(81, 79)(69,67)(81,69)(79,67) 




(21, 16)(18, 13) 


160 


(88, 84)(76, 72)(88, 76)(84, 72) 


36 ■ 


(21,21)(15,15)(21,15) 


*162 


(90,81)(81,72) 


40 


(24, 22)(18, 16) 
(24, 18)(22, 16) 


1 64 


(91,83)(81,73)(91,81)(83, 73) 


*50 


(30, 25)(25,20)(29,28) 
(29, 22)(28,21)(22,21) 


*1 70 


(94, 87)(83,76)(94,83)(87,76) 
(92, 91)(79,78)(92,79)(91,78) 


52 


(31,27)(31,25) 

(25,21)(27,21) 


*178 


(97,94)(97,84)(84,81)(94,81) 


*58 


(34,31)(27,24) 

(34,27)(31,24) 


1 80 


(99, 93)(87,81)(99,87)(93,81) 


64 


(36, 28)(36,36)(28,28) 


* 1 94 


(106, 101)(93, 88) 


68 


(39,37)(31,29) 




(106,93)(101,88) 




(39. 31)(37,29) 


196 


(105, 105)(91,91)(105,91) 


72 


(42,36)(36,30) 


200 


(110, 100)(108, 106)(100,90) 


*74 


(43,38)(36,31) 
(43, 36)(38,31) 




(94,92)(94,108)(106,92) 


80 


(46,42)(38,34) 


All pos sible 


complementary 


*82 


(46, 38)(42,34) 
(46,45)(37,36) 


code lengths up to 200 with 
unordered possible weiehts. 








(46,37)(45,36) 


*90 


(51,48)((51,42) 
(42, 39)(48,39) 




TABLE 2„1 


00 


(56,49)(49,42) 





100 (55, 55)(55,45)(45,45) 

(57, 51)(57,49)(51,43) 
(49,43) 
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The shortest possible complementary pair is A=ll, B = 10, This pair or 
any of its transformations is called a kernel of length two, or the quad. 

A kernel is a basic length code which cannot be decomposed into shorter 
length codes by an inversion of the standard generating methods to be 
explained later. Some possible kernel lengths are 2, 10, 18, 26, 34, 

50, etc. although codes for all of these do not exist. Complementary 
sequences which are not kernels are called composite complementary 
sequences. 

The above list of kernel lengths did not include n=4 or n = 8 
which are possible complementary sequence lengths. These if they 
exist must therefore be composite. All possible codes of length four 
can be generated from the exhaustive list of possibilities for four binary 
digits given in Table 2. 2. 



0000 


1000 


0001 


1001 


0010 


1010 


0011 


101 1 


0100 


1 100 


0101 


1 101 


0110 


1110 


0111 


nil 



Table 2.2. All Possible Binary Numbe rs of Length Four. 

The unordered possible pairs of ones in the codes are (1, 1), (3,3) 
and (3, 1). Thus all possible codes of length four have either one 1 or 
three Ps because of this limitation. This reduces the table to the 



eight following numbers of length four: 



0001 






0111 


0010 






1 101 


0100 






101 1 


1000 






1110 


The possible unordered (1, 1) 


code 


pairs are the following four: 


1000 


1000 


0001 


0001 


0100 


0010 


0100 


0010 
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The possible unordered (3,1) pairs of codes are the following eight: 

0111 0111 1110 1110 1011 1011 1101 1101 

0010 0100 0100 0010 0001 1000 0001 1000 

The possible unordered (3,3) code pairs are 

0111 0111 1110 1110 

1011 1101 1011 1101 

Comparing the last code to the kernel of length two, A=1 1 , B=10, 
shows that it might have been constructed by writing in time sequence 
S = AB=aj^a2bj^b2 = 1110, 

= AB=a^a2b^b2 = 1101. 

The next to last sequence pair might have been formed by inter- 
lacing A and B in the following manner; 

Ti — aj^bj^a^b^ — 1110, 

T 2 = a^Bj^a^^b^ = 1011. 

All the rest of the 16 pairs could be considered either as transforma- 
tions of the last pair, or as being formed by the same two composite 
generating operations used above on the 8 transformations of the ker- 
nel of length two. A general proof will now be given to show that both 
the time sequence scheme (Si, S 2 ) and the interlace scheme (Ti, Tg) 
will always form complementary sequence pairs providing that A and 
B are a complementary sequence pair. 

Theorem 2o 7 

If A.— a-ia^^a^. ...... a 

l c 3 n- 1 n 

•— b-ib-^bQ. . .. .. ©b 

icy n- 1 n 

are a complementary sequence pair, then 
C — ■ a.a'^.c.a bib'^....... b 

1 Z n 1 Z n 

** ^n^l ^Z* * * ® ® ^n 

are a complementary sequence pair. 
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Theorem 2, 8 



If (A, B) are a complementary sequence pair, then 

*“ a 1 b 1 a -jb eo« a b 

i i Z Z n n 

D = 

form a complementary sequence pair. 

The following notation will be used in the proof: 
i=n- j 

U = S a.® a. . for all j 1 < j < n-1 

^ i “ 1 ^ J 

L = ® 1 for all j a < j < n-1 « 

a El i+j J — J - 

i=l 



It is understood that the two CD sequences in the two theorems 
above are not the same but because the proofs are identical two sets of 
symbols will not be used. For.^each spacing j there are three possible 
ways for bits to match, within A, within B, or from a bit of A to a 
bit in B. These possibilities are denoted by and 

U^, respectively*, 

The necessary and sufficient condition for C and D to be a 

complementary sequence pair is given by 

= U +Ur + U - where L U .. 
d a b ab c d 

Now L = U, and L, = U , i r- ... r i 

a b b a by deiinition 01 complementary. 



Adding 



+ hi = 



but 

Ut_ = Ur since 
b b 

so substituting 



i=n-l i=n-l 

E b^©b . = E 

i=l 1=1 



L + JU = U + Ur 

a D a b 

All bit matches from A to B in the C sequence are paired in the D 
sequence with bit matches from A to B therefore. 
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L U = U r , 
ab ab 



and adding 



L + L,+L,=U +Ur + Ur 

a D ab a b ab 



or L = U , 
c d 



Golay describes two other methods of generating composite sequen- 
ces from shorter complementary sequences. These methods will not be 
used too extensively in the rest of the paper and will therefore just be 
mentioned rather than proved. Given two complementary sequences 
pairs (A, B) (C,D), A of length m and C of length n. 

Let 



Ui = A^^A^? A^n B^^ 

U3=A^’^.... ....A^^ “B^n 



(2.6) 



and 

Vi = A*^^ B'^^ A*^n B^n 

a - j - (2. 7) 

V 3 = A n B^n A ^ B^\ 

It can be shown that the pairs of U are complementary as are the pairs 
of V. Where, if an exponent is a one the A or B code is left un- 
changed, and if the exponent is a zero the A or B code is comple- 
mented. The lengths of codes which can be generated from these two 
methods are 2mn where m is the length of the (A,B) sequence pair 
and n is the length of the (C , D) sequence pair. There is also another 
special method of generation which applies only to codes of length 2 . 
This method is explained in Chapter 5 where it is necessary for the 
completeness of a proof. 

The kernel of length 2 , or the quad, has been discussed in the 
previous pages. There are two kernels of length 10. These are 
1001010001 and 1000000110 , and a second kernel 0101000011 and 
00001 001 10. The next possible kernel size is 18 and Golay has proved 
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by exhaustive search that no kernel of this length exists^ Kruskal has 

1 3 

since then completed the proof analytically. One of the goals of the 
writer in his study was to make an exhaustive search for codes of kernel 
length Z6, This search disclosed there was only one kernel of length 
26, Chapter 7 of this paper describes this search in detail and 
Chapter 8 describes attempts to find kernels of length 34. 
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CHAPTER UI 

THE GROUP OF OPERATIONS 

A large portion of this investigation is devoted to the set of allow- 
able transformations on complementary sequences. These allowable 
transformations possess the characteristic that the length of the sequen- 
ces and the complementary property are invariant under transformation. 
To be useful in searching for new kernels it is also desirable that these 
transformations have the following features: 

1. That they form a group for all general transformations. 
Where a general transformation is one which can be 
applied to any complementary code pair. 

2. The transformations as represented be simple manipula- 
tions on the elements of the group. 

3. The supplementary characteristic (to be explained in 
Chapter 6) be a property of the representation. 

4. The parity check be easily made. 

Some of these transformation operations were discussed, and the 
proof of invariance given for all group generators, in the background 
chapter. Chapter 2, on complementary sequences. As a continuation, 
a brief review of the transformation operations and a tagging of symbols 
to these operations will prove useful. Each basic operation symbol and 
its definition is given in Table 3. U More definitions are given than are 
necessary to generate the group; however, the redundance seems to the 
writer to be an aid to understanding. 

A check through Table 3. 1 shows that the transformations affect 

the complementary sequence pair in four different ways. Ci , Cg , 

Ti , Ts affect all the bits of one of the code pair. affect all bits 

of both codes, while Ai and Ag affect half of the bits of both codes. The 

identity, I, of course has no affect on either code. A reasonable way 

to form the elements of the group to account for the above affects would 

be to divide each code of the pair into two units. These would be A^ 

(even bits) and A^ (odd bits), also (even bits) and (odd bits). 
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The Identity 



a. a. 
1 1 

b. -* b. 
1 1 



A code time reversal 



a. -a ,, . 

1 n+l-i 

b. -b. 

1 1 

B code time reversal 



a. 

1 



a. 

1 



i n+l-i 

Time reversal of both codes 



a. a , , . 
1 n+ 1 - 1 



Ai 



Ag 





^2i ^ ^2i 


^2i-l 


^2i “* ^Zi 



b. - b . . ; 

1 n+ 1 - 1 

Alte r odd bits 
a 
b 

Alte r even bits 
a 
b 

Complement the A code 

a. “* a. 

1 1 

b. - b. 

1 1 

Complement the B code 



2i-l*' ^2i-l 


^ZC ^2i 


2i- 1 ^2i- 1 


b, b,. 

2i 2i 



a. 

1 



a. 
1 

b. 
1 



Complement both codes 



a. 

1 



a. 

1 



Some 

Exchange the A code with the B code 



b. 

1 



b. 

1 



a. 

1 



b. 

1 

a. 

1 



operations. 
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TABLE 3. 1 




i 



y 



■A. — 1 

o 13 d n- 1 

A. = 3,^3, .3./ , *.o.o«.,3. 

e Z 4 0 n 

1 

o 13 d n- 1 

e Z 4 o n 

However, this particular division for representing the complemen- 
tary sequences proves to be very awkward for the time reversal opera- 
tions (Ti ,Ts , T) and for the ease of the parity check condition. 

Fortunately, one simple change, that of time reversing the second and 

14 

fourth row (even bits of A and B), clears up these difficulties. 

The method for converting from a complementary sequence pair 

into the standard group form is therefore 

I— a^a^ar:...«e3< i 

13 5 n- 1 

-^2 (3.1) 

III=b.b,b(- b , 

B=b,b,b. b ,b ( 13 5 n- 1 

1 ^ ^ n- 1 n \ j.r. . . 

^^=Vn-2*-”--^2 

H the matrix (I II III IV) is taken as the symbol for the above grouping 
the Identity operation, I, should yield (I II III IV). Post multiplying 
(I II III IV) 

10 0 0 

by 0 0 10 (1 11 111 1^)’ therefore I can be taken 

0 0 0 1 

as the identity matrix of rank 4, using the operation matrix multipli- 
cation on the right. Similarly Ti has the effect of exchanging the first 
and second columns yielding (II I III IV), Now post multiplying 



(I II III IV) by 

0 10 0 , 

0 0 10 yl®l<3s (II I III IV) the desired result so the 

0 0 0 1 

matrix just used could represent Ti . Similarly, 
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yi IL 



T= 



10 0 0 
T _ 0 1 0 0 

0 0 0 1 

0 0 10 



0 1 0 
1 0 0 
0 0 0 
0 0 1 



Ai- 



1 0 0 0 1 0 
0 1 0 0 0 ! 
0 0 i 0 ^2 Q Q 

0 0 0 1 0 0 



0 1 0 0 0 1 

0 r - - 0 1 0 0^_0 

1 0 0 1 0 0 

0 0 0 0 1 0 



0 0 0 0 1 0 

0 0 0 0 0 1 

1 0 ^ " 1 0 0 0 

0 i 0 10 0 



0 0 0 
1 0 0 
0 1 0 
0 0 I 



Matricies are often a very convenient method of expressing trans- 
formations. However, when the elements are few and the numbers in 
them simple it is often more convenient to express the operation dif- 
ferently. For instance (I II III IV) operated upon by Ti yields 

(II I III IV). This could be just as well expressed as 

I II III IV 
^ " II I III IV • 

This indicates that the top row is transformed into the bottom row by the 
operation, or in even more shortened form, Ti = (II I III IV) with the 
top row understood which means, as is common practice, Ti (I II III IV) 
= (II I III IV). Another example is Ai = (I II III IV), [which means 
Ai (I II III III) = (I II III IV)] where the bar across I and III indicates 
the negative element in the matrix or in actuallity the complement of the 
portion of the code contained in I and III. Table 3, 2 is a list of all 
32 transforms generated from Table 3. 1 both by symbol and group 
element representation. There are actually 64 transform operations 
in the entire group. However, the exchange operation (III IV I II) is 
not considered along with all the 32 elements it would generate. In 
this paper the pair (I II) and (III IV) are generally considered as 
being an unordered pair, which eliminates the exchange operation. 

The second feature listed under desirable characteristics, that 

of simple manipulation of the elements of the group is satisfied. 

The supplementary property which is discussed in Chapter 6 is 

also satisfied by the formulation of groups from complementary 

34 



sequences by this method. The parity check, a ®b.®a .®b . = 1, 

’ 11 n+l-i n+l-i 

is also satisfied under all transforms and is easily checked since each 
column of equation 3. 1 is one solution of the parity check. 



I = a 
II = a 


1^3 


^5 * * 


e • <» ^ • • • e • • 

1 n- 


■ 1 










3i 

n n 


-2 • • 


‘ • •‘^n+l-i ^2 












III = b 


1^3 


b ^ . 


*««b. «••• b 
1 n- 


1 










IV = b 


b 

n n 


-2 • • 


• ” ''^n+i-*l’^2 












I = ( I II 


III 


IV) 




R 




(I TT 


IV 


m 


Ai= (T II 


Hi 


IV) 




Q 


= 


(u T 


IV 


ST) 


A2= (I TT 


III 


iv) 




P 


= 


(IT I 


Tv 


UI) 


Ti= (II I 


III 


IV) 




O 


= 


(n r 


m 


IV) 


Ts= (I II 


IV 


III) 




N 


= 


(U I 


m 


IV) 


T = (II I 


IV 


III) 




M 


= 


(T n 


IV 


III) 


Ci= (T II 


III 


IV) 


1 


L 




(in 


m 


W) 


C3= (I II 


m 


W) 




K 


= 


(1 Ti 


Tv 


m) 


C = (1 TI 


m 


W) 




J 


= 


(IT 1 


Tv 


m) 


z = (II T 


TTi 


IV) 




H 


= 


(I II 


Tv 


ITir 


Y = (IT I 


m 


IV) 




G 


= 


(ITT 


IV 


III) 


X = (IT I 


III 


IV) 




F 


= 


(I "n 


iiT 


IV) 


w = (II T 


III 


Tv) 




D 




(I II 


m 


Tv) 


V = (1 II 


IV 


iTi) 




0 


= 


(II I 


Tv 


ITT) 


U = (I II 


fv 


III) 




B 


= 


(IT I 


IV 


In) 


s = (I H 


TT 


III) 




IT 


= 


(II T' 


IV 


lU) 



E = (III IV I II) 

Elements of the unordered operations group. 



TABLE 3. 2 
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TABLE 3.3 

Operations Multiplication Table 
(Sicle) (Top) = Operation 



Another feature inherent in this formulation of the group, which is 
very useful in searching for new kernels, is the interlace breakdown 
into the group. This is useful because a kernel may be considered as the 
interlace of two half length codes where, as will be shown in Chapter 7, 
the number of zeros and ones in each half length code are easily pre- 
dete rmined. 

Although this group of operations is indeed a group, it has never 

been proved formally. In order for a collection of elements to be a 

1 5 

group, they must satisfy the following four conditions: 

1. A group has closure. 

2. A group is associative under the group operation. 

3. A group has an identity element. 

4. Each element in a^ group has an inverse. 

That this operation group has closure is demonstrated by Table 3. 3 
which is the full multiplication table of this group of order 32. The 
elements under the operation multiplication are associative since they 
are expressible as matricies and matrix multiplication is associative. 
The element I is the identity, and as seen from Table 3.3 all elements 
have an inverse. Therefore these 32 elements have satisfied all the 
requirements for a group. 

There are three special subgroups of order 4 which within the 
subgroup deal with only one type of operation, Ti Tg TI, Ai As Cl, 

Cl Cs Cl. All of these subgroups have the unusual property that each 
element in the subgroup is its own inverse. This group of order 32 
has been identified as isomorphic with Senior’s group number 44. 

This will be discussed in some detail at the end of this chapter. 

A study of complementary codes in (I II III IV) form, or as we 

will call it sequence quadruple form, reveals many interesting proper- 
ties. Some of these properties and their formal proofs will follow. 
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In these proofs the complementary pair (Aj,B) will be assumed to be of 

i=j 

length n. The symbol L,^ will mean S ‘ ® 

j. 1 < j < n-1 and similarly Ug = ^ 

i=l i ^+i-j * 

Theorem 3. I 

If (AjB) are a complementary sequence pair and if either (I, II) 

or (III, IV) in sequence quadruple form are complementary, then the 

other pair must also be complementar y» 

1. Assume (I, II) are comlementary„ 

Lij = Ujj by the definition of complementarity. 



2 . 

3. 

4. 



Uj = Ljj by the definition of complementarity. 

Adding the equations in steps 2 and 3 gives 

L, + = U, + U„. 

Since only even spacing s are concerned in the sequence 
quadruple form, when related back to the original (A,B) 
pair, the likes of A are equal to the unlikes of A for all 
even spacingSo This will be expressed as for 

even space. 

Li^=Ug and by the definition of complementarity. 

for even spacings. 

Jo Jj 



is equal to the likes plus the unlikes, and is dependent upon j. 



6. 


La=Ub and 


7. 


Therefore ] 


8. 






Ill III J 




is equal to t 


9. 




1 o„ 


Adding 8 a 




(UlII+Uiv) = 


1 1. 


bii+^iv ' ^ 


1 2. 


hv-Uiir 0 




in step 1 1. 


1 3. 


Therefore 



7 -^Iir^IIl' 

'j “ 






Lijy ”^iJi (Iffs ^V) pair is complementary. 

Since the proof would have been the same if (Illg IV) were 
chosen complementary rather that (I, II), the general state- 
ment is true. 
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Theorem 3o 2 



If (A,B) are a complementary sequence pair in sequence quadruple 
form and if one of the pair (I, II) or (III, IV) is not complementary 
neither is the other. 

1. Assume (I, II) not complementary. 

2. There are two possibilities for (III, IV) either complementary 
or not complementary. 

3. Assume (III, IV) are complementary, then by the Theorem 3.1 
the (I, II) pair must also be complementary, but this is con- 
trary to assumption. 

4. Therefore (III, IV) cannot be a complementary pair. 

>;c 



Theorem 3. 3 

A necessary and sufficient condition for (I, II) and (III, IV) to be 
two complementary pairs when a complementary sequence pair is 
written in the sequence quadruple form is if the likes for even spacing of 
one of the code pair is equal to the unlikes at the same spacing for the 
same code of the pair. First the sufficient portion: 



1 . 

2 . 

3, 

4. 



6 . 

7, 

8 . 
9. 



L, = U« for even spacings. 

A A 

Ll + Lj-j. = L.^ for even spacings. 

Uf + Ujj = - L.^ for even spacings. 

L.J + Lj-j + Uj + Ujj = even spacings by adding the 

equation in step 2 and step 3. 

L<i + Uj = kj since the total of likes and unlikes must equal the 
pos sibilities. 

Ljl + Ujj = kj . 

Lij + Uj + Ltjj 4 Ujj = 2k. = by adding the equations in steps 

5 and 6 and comp>aring to those in step 4. 

Lj + J_jj = kj rewriting step 2 with k^ substituted for 

Ljj - Uj = 0 subtracting the equation in step 5 from that in 
step 8. 
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10. Therefore hi = Uj and (I, II) are a complementary pair and by 
Theorem 3. I (III, IV) are also a complementary pair. 



11 . 

12 . 

13 . 



For the proof of the necessary condition assume that (I, II) are 
complementary and that for some even spacing. 

Therefore Uj + for some even spacing. 



However, by a ssumption (I, II) 

Lj = Ujj and Uj = Ljj. 



are complementary and 



14. Adding the two equations in step 13 gives . 

15. This is incompatible with step 12, therefore it is necessary that 

Ua = or Ug = Lg for all even spacings in order to have 

(I, II) and (HI, IV) complementary. 












Theorem 3.4 

If the sequence quadruple form of a complementary pair can be 
transformed into the pattern (I II I TT), either the original pair is the 
quad or (I, II) and (II, IV) are complementary pairs, 

1. Given that (I II HI IV) is transformable into (I II I II), and 
is so transformed. 

2. I-'i‘*’I-'ii‘*'I-'in^^Iv”^I^^II^^III^^IV~ supplementary 

property, Theorem 6.1. 

3. I-'ii=I=ii' ^ir^TT ^i®^i+j“ ^i®^i+j Theorem 2. 1. 

4. Therefore 2L,j+2I_,jj= 2Uj+2Ujj= 2k^ by substituting into step 2, 
and using the condition of step 1. 

5. k. . 

6 L,.+U,= k. since the total of likes and unlikes must equal all 

I I J 

pos sibilitie s. 

7. I-'ii"I-^i= ® subtracting the equation in step 6 from that in step 5. 

8. Therefore Lg= Uj and (I, II) are complementary. 

9. Since (I, II) are complementary, so are (IjlT) and any general 
transformations of (I, H) and (III, IV). 



This last property should be extremely useful in decomposing long 
codes, since it avoids the tedious task of checking to see if (I, II) are 
a complementary pair^ Whenever a complementary pair expressed in 
sequence quadruple form has the pattern (I II I TT) , certain additional 
transformations are allowed which are not true in general These 
allowable operations are the time reversal of I and III, II and IV, 
or both. Also allowed as an operation is the altering of either even bits 
or odd bits for all four sequences. The proof of all these operations 
being allowable is quite simple, but does require an addition to the no- 
tation to cover these particular operations. Double or quadruple sub- 
scripts will be used to denote these operations, T 2^^ and ^^^234 

represent respectively the time reversal of (I, III), {II, IV) and the 
time reversal of both pairs. represents the altering of the odd 

bits in all four sequences while stands for the altering of even bits 

in all four sequences. Actually complementing any two, or all four of 
the four sequences will also leave a complementary pair when reassem- 
bled, but each of these operations is already accomplished by the 
general transformations. 

Theorem 3. 5 

Whenever a complementary sequence pair is put into sequence 
guadruple and then transformed so that the resultant group pattern is 
(I II I TT), I and III, or II and IV, or both, can be time reversed 
and the sequence pairs formed on reassembly from thequadruple, after the 
operations or will still be complementary sequence 

pairs. 

1. Given, (A,B) is a complementary pair of length n transform- 
able in quadruple form to (I II I IT), and is so transformed. 

2. (I, II) is a complementary sequence pair by Theorem 3.4. 
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I 



'=cS9i*» 



I 



' 



3. Therefore (1^,11), {^H) are complementary pairs 

by Theorem 2. 3, 

4. From this it is seen that the original sequence quadruple after 

the operations ^1234 ^ J 

1} ^ 11) and (J_ ^ i. ZD* ^hen interlaced in the customary 
manner to form pairs, the pairs they form will be complementary 
sequence pairs by Theorem 2« 8 ^ 

>;< >;< ^ 



Theorem 3, 6 

Whenever a complementary sequence pair is put into sequence quad- 
ruple and then transformed so that the resultant group pattern is (1 11 
I II), the operation the altering of odd bits, or the operation 

A 241 the altering of even bits, will allow the quadruple upon reassembly 
to still be a complementary pair. 

1. By hypothesis A and B are a complementary pair which 
when expressed in sequence quadruple form are transformable 
to (I II I TT), and is so transformed. 

2 . (I, II) is a complementary sequence pair by Theorem 3.4. 

3. (I, II) when operated upon by either or A^ is still a 

complementary pair by Theorem 2.4„ 

4. If the symbol ( 1, 2) is used for the pair (I, II) altered it is 
seen that the pair formed from the interlace of (12 12 ) is 
a complementary pair by Theorem 2. 8 . 

5|C 5jc 

These last two theorems should prove quite useful when trying to 
estimate an upper limit on the possible number of composite comple- 
mentary sequences of any particular length. An example to illustrate 
the use of Theorem 3.5 might prove helpfuh Taking pair number one 
from Appendix II, 
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3 




‘I ii II II I 



A=11010001 1101 1110 
6=1000010010001011 
1 = 10001011=1 
11=01 111011 = 11 
111 = 10001011=1 
IV = 1 0000100 =ir 

Since this is already in (I II I II) form without transforming it, 

can be applied directly giving 

1=1 1010001 
11=01 11 1011 
111=11010001 
IV = 10000100. 

Reassembling gives 

A^=1111001101010110 

6^=1010011000000011 

which is identical with number three in Appendix II. 

As an example, to illustrate the use of Theorem 3.6, applying 

to the same (A, 6) pair gives 

1. =11011110 

2. =00101110 

3. =11011110 
I . = 11010001 , 

reassembling give 

A^=101101 11 10111000 

6 ^= 1110001011101101 . 

This is the same as number 6 in Appendix II with the transformation 
TC 2=0 applied. 
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The (I II I II) form, or one of its transforms, appears quite 
often in the study of composite complementary sequences. Since two 
of the components are the same, and the other two are complements, 
it is reasonable to assume that some of the transforms are redundant, 
that is yield indistinct results. Table 3.4 lists all the transforms 
with the elements 1, II and their complements. If (A, B) are taken 
as an ordered pair, the number of different codes formed from the 
transformation group are 3Z instead of 64 because A 2 (I II I II) 
and E(I II I Hy are the same. If (A,B) are taken as an unordered 
pair the number of transforms forming different codes is decreased 
from 32 to 16. Table 3.4 lists each transform with its unordered 
mate. 

The only possibility of having less than 16 unordered pairs of 
codes which are distinct would be if the group form could be transformed 
into the pattern (III I). The only case known to the writer where 
this is true is the quad. For one version of the quad this is 

A=ll B=10 

1=1 

11=1 

III=1 

IV=0. 

Table 3. 5 shows the eight ordered pairs, or if the cross index column 
is used the four unordered pairs, for the quad. The four possible code 
transformations for the quad are: 

12 3 4 

A= 11 11 00 00 

B= 10 01 10 01. 



45 



Patte rn 



Transforms with pattern Cross index number 
for ordered pairs for unordered pairs 
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T) 


Ci.O.K, J 
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I) 


C, M, L,G 
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(I 
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Z,R,Q,F 


f. 


8. 


(I 


I 


I 


I) 


X,U,P,D 


4. 



Redundance due to (I I I T) . 
TABLE 3.5 



This exhausts the transformations of special types and the redun- 
dancies caused by (I II I II), The operations group will now be 
compared to Senior’s group number 44 of length 32 as was mentioned 
earlier in this chapter. 

The identification of any group with a group in Senior’s list, is 

accomplished through establishing a one to one correspondence between 

16 As will be demonstrated, Senior’s group num- 
group generators, ^ ^ 

ber 44 of order 32 is in one to one correspondence with the unordered 
operations group of complementary sequences expressed in sequence 
quadruple form. Therefore the two groups are isomorphic. 

Senior’s generators will be given in lower case letters, while those 
of the operations group will be in their standard upper case symbols. 
Five generators are required for Senior’s number 44; let these be 
a,b,c,d,e. The following relationships between generators are re- 
quired for Senior’s 44: 

a^=b^=c^=d^=e^ = 1 where 1 is the identity, 

b"kb=a, c"^ac=a, d”^ad=a where x ^ is the inverse of x, 
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1 . 

2 , 



W‘ 




1 



3 . 


c ^bc=b. 


d ^bd=b. 


4. 


d ^cd=c. 




5. 


e ^ae=a. 


e ^be=b, 


6. 


e ^ce=ca 


, e Me=bd 



A check of the above requirements shows that a,b,c,d form an 
abelian group of order 16 while the ^’e'* generator is commutative 
only with a,b and is non- commuta tive with c,de 

Let a=Ci, b=C 2 > c=Ti, d=T 2 , e=Ax » From Table 3*3 it is seen 
that the following relations are true* 

1. C? = C|= T? = T|= A?=I 

2 . C gC iC2 = Cl, TiCiTi = Cl , T2 CiT 2 = Ci . 

3. Ti C2 Ti = C2 • 

4. T2 Ti T2 = Ti * 

5. AiCi Ai = Cl , Ai C 2 Ai = Cs * 

6. Ai Ti Ai = Ti Cl = 0, Ai T 2 Ai = C 2 T 2 = H . 

This one to one correspondence between the generators of Senior^s 
group number 44 and a set of generators from the unordered operations 
group on complementary sequences is necessary and sufficient to prove 
the two groups are isomorphic* 

The diagrams for Senior's group number 44 of order 32 and 
those for its subgroups were used as a guide line to check the generation 
of subgroups from the operations group of unordered complementary 
sequence pairs. These subgroups of order 1,2, 4, 8 and 16 are shown 
in Appendix IV. All groups that were shown in Senior's diagrams are 
included in these charts; however. Senior's diagrams included only 
normal subgroups and, although much care was taken in the generations 
of the subgroups of the operations group, there is a possibility that 
Appendix IV is not exhaustive due to some undetected non- normal sub- 



groups. 
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Senior in conjunction with others has rewritten his work and 

17 

changed his designating method. Under this new method his 

number 44 of length 32 is now The ordered operations group 

of length 64 is also identified as 

An identification of the group with Miller’s number 23 of order 

32 was also accomplished for those who prefer the substitution group 

1 8 

method of identification. 
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CHAPTER IV 



HAMMING DISTANCES OF COMPLEMENTARY SEQUENCES 

The Hamming distance is a well known property of certain classes 
of codes. This very useful property has had most of its applica- 

tions in the fjelds of binary error detecting and error correcting codes. 
This chapter will demonstrate the usefulness of Hamming distances in 
the field of complementary sequences, as both the code pairs and their 
sequence quadruple formulation have certain invariances in Hamming 
distance under various conditions. 

The Hamming distance, D(U, V), of two binary vectors or sequences 
U and V is defined as the number of positions in which these'two binary 
vectors differ. 

i=n 

D{U, V) = E u. ® V. (4. 1) 

i=l " 

For example, the Hamming distance of the two binary vectors 

U=1011001 

V=0111001 

7 

is D(U,V) = E u. ® V. = (1 + 1+O+O+O+OH-O) = 2 
i=l 

The Hamming weight or just weight of a binary vector is its distance 
from the null vector. 

Theorem 4. 1 

All complementary sequence pairs of length n have a Hamming 

distance of n/ 2 . 

By hypothesis 

TV— a 1 a -5 a ^ . . « « « • ^ 

12 3 n 

and 

B ““b-b-jbo* o «• e • b 
12 3 n 

A and B form a complementary sequence pair. 
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r:'SSS(S*S* 






1 . 



2 . 



3. 

4. 



All complementary sequence pairs satisfy the parity check 



a.«a 
1 n-i+1 



= 1 , f or all 



i n-i+1 
i, 1 < i < n/2 by equation 2, 4» 

Since mod 2 addition is commutative 

a. ®b. ®a . ®b = 1 . 

1 1 n-i+1 n-i+1 • 

Therefore either (a. ®b.) = I or (a ®b .,i), but not 

11 ' n-i+1 n-i+1' 

both. 

For each i, 1 < i < n/2 one of the two vector positions in 
step 3 will have a one and the other vector position a zero. 

Summing up over all possible i, gives n/2 for the 
Hamming distance. 



^ jje 

When a complementary sequence pair is decomposed into a sequence 
quadruple, each pair of sequences in the quadruple, which form one of 
the complementary pairs, has a Hamming distance which will be shown 
to be characteristic of the code length. The sum of the two pairs of 
Hamming distances from the quadruple is an invariant and is equal to 
n/2 the Hamming distance of the complementary sequence pair. 

For example. 



A=1001010001 
B=10000001 10 



are a complementary sequence pair of length 10. 

Decomposing them into the standard sequence quadruple, (1 11 111 IV), 
give s 



I 

II 

III 

IV 



10000 

10110 



= 10001 
= 01000 
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I ® II =: 00110 D(I, U) = 2 

III ® IV = 1 1001 D(III, IV) = 3 

I ® II © III ® IV = 11 111 D{I«n, TII®IV) = 5 . 

This distance is half the length of the complementary sequence pair, 

and is required because of the necessary condition of the parity check. 

If the distances (2,3) for the pairs given in the above example are 

taken as an unordered rather than an ordered pair, this (2,3) pair 

will be invariant through all general transformation operations. 

All of the general transformations on the sequence quadruples 

have one or more of the following properties: 

1. Change the order of (I, II) or (III, IV) or both, 
for example, Ti , Ts, T . 

2. Complement the pairs (I, II) or III, IV) or both, 
for example, Ci ,Cg , C . 

3. Complement one from the pair (I, II) the other from 
the pair (III, IV), for example, Ai , A^ „ 

4. Exchange (I, II) for (111,1V), for example, E. 

Since all operations in the group can be generated by multiplication 
from the set listed above, an examination of the effect on Hamming 
distance by these transformations will be adequate to prove the invar- 
iance. 

Theorem 4. 2 

A change in order of (I, II) or (III, IV) or both will not change 

the Hamming distances of the pairs (I,II)^ or (III, IV)« 

n 

1. By definition D(U, V) = S u^ ® v^ 

i“ 1 

2. u. ©V. = V. ®u. modulo 2 addition is commutative 

11 11 

i=n i=n 

3. D(U,V) = S u ©V. = Z V. ®u =D(V,U) 

i=l ' ' i=l ' ' 



51 



Theorem 4. 3 



Complementing (I, II) or (III, IV) or both will not change the 
Hamming distances of the pairs (I, II) or (III, IV). 

1. u. ® V. = u. ® V. by Theorem Z. I 

2. Therefore S u. ® v. = Z u. ® v. 

11 11 

3. and D{U,V) = D(U, V) 

5 {< 

If it is assumed that the length of the sequence quadruples formed 
from the complementary sequence pairs have a length K and that the 
distance of the first pair is R, then the following theorem can-be proved: 

Theorem 4. 4 

Complementing one of the pair (I, II) and one of the pair (III, IV) 
results in the exchange of their Hamming distances, or expressed dif- 
ferently, gives the same distance pair in opposite order. 

1. Let D{I,II) = R, therefore 

D{III, IV) = K-R 

as the parity check requires D(I®-II, III'®IV)=K 

2. Complement one of the first pair in the (I II III IV) form, 
and also complement one of the second pair. 

Each vector position that was formerly alike and summed 
to zero is now different and sums to one. 

Similarly each vector position that was formerly different 
and summed to one is now alike and sums to zero. 

3. D{T, II) = D{I,II) = K-R from 2. 

4. D{in, IV) = D{III,TV) = K-(K-R) =R from 2. 

>!« ^ 

Theorem 4. 5 

The exchange operation has the effect of changing the order of 
Hamming distances of the pairs (I, II) and (III, IV) when expressed 
in sequence quadruple form. 



52 



1. The pairs (1,11) and (III, IV) are not changed in any way 
except by position in the exchange operation. 

2. The only way to effect the Hamming distances of a pair of 
vectors is to change some of the components of the vectors. 

3. Therefore the two weights are the same but reversed in 
order. 



An examination of the Hamming distances for the known kernels of 
2, 10 and 26 is now in order as these distances are a characteristic 
of the kernel, as are the vectors formed in the process of obtaining the 
distances by the modulo 2 sums. 

The first case to examine is the trivial case for the kernel of 
length 2. 

A=ll B=10 

Broken down into sequence quadruple representation 

1=1 D(I,II)=0 

11=1 

III=1 D(III, IV) = 1 

IV=0 

There are two kernels for length 10 to be examined. 

A=1001010001 
B=10000001 10 



1=10000 


1 ® 


11=001 10 


D{1, 11)= 2 


11=101 10 








111=10001 


111 ® 


1V=11001 


D(lll, 1V)= 3 


1V=01000 








A=0101000011 
B=00001001 10 








1=00001 

11=10011 


1 « 


11=10010 


D(l, 11)= 2 


111=00101 


111 § 


1V=01 101 


D(lll, 1V)= 3 
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Both kernels of length 10 have the unordered Hamming distances of 
(2,3) when expressed in sequence quadruple representation^ However, 
their ’’Hamming vectors", the unordered pairs (001 10,1 1001) and 
(10010,01101) are different and are characteristic of the kernel of 
their origin under all possible transformation operations. Hamming 
vectors will be considered at length in the next chapter and proof of the 
above statement will be deferred until then. 

There is only one kernel of length 26 (Chapter 7) and it is 

A=01001 101 1110101 11100111010 
B=1011 00 100001 1111 1100111010 

1=0010111110111 I ® 11=0000001101010 D(I,II)= 4 

11=00101 10011 101 

111=1101001110111 III ® IV=1111110010101 D(III,IV)= 9 

IV=0010111100010 

This is a complete listing of the Hamming distances for all known 
kernels and it is exhaustive for n=2,10,18 and 26. It is worthwhile 
to note that the sum of the unordered pairs of distances is always equal 
to one half the length of the code and that the Hamming vectors of the 
pairs are complements. Both of these conditions are due to the parity 
check being a necessary condition. 

A natural extension of the Hamming distances of kernels is the 
consideration of Hamming distances of composite codes when put into 
sequence quadruple form. Six different code pairs of length eight are 
listed in Table 4. 1. These pairs will be examined for a possible in- 
variance. After this invariance is noted a general proof will be given 
for the characteristic distance of composite codes. 

The first four examples in the table when decomposed into sequence 

quadruple representation have (I, II) and (III, IV) as code pairs which 

are still complementary. These codes are half the length of the original 
pair, and being complementary, have a Hamming distance one half of 
their length, or n/4 as compared to the original codes of length eight. 
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(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



A=11111001 

B=11001010 

1=1110 

11=1011 

111=101 1 
1V=0001 

A=ll 110110 
B=11000101 

1= 1 1 0 1 
11=0111 

111=1000 

1V=1101 

A=11101101 

B=11100010 

1=1110 
11= 1 1 0 1 

111=1 101 
1V=0001 

A=1101 1110 
B=11010001 

1=1011 

11=0111 

111=1000 

1V=1011 

A=11101011 
B=ll 100100 

1= 1 1 1 1 
11=1011 

111=1100 

1V=0101 

A=1011 1110 
B=10110001 

1= 1 1 1 1 
11=01 10 

111=1100 

1V=1010 



H=(00110011) 
1€ 11=0101 

111 ® 1V=1010 

H=(00110011) 

1 ® 11=1010 

111 ® 1V=0101 

H=(000011 11) 

1 ® 11 =001 1 

111 ® IV = 1 100 

H=(000011 11) 

1 ® 11 = 1 100 

111 ® IV = 0011 

H=(000011 11) 

1 ® 11 = 0110 

111 ® IV = 1001 

H=(000011 11) 

1 ® 11=1001 

111 ® 1V=0110 



D(l, 11) = 2 
D(lll, IV) = 2 



D(l, 11) = 2 
D(111,1V) = .2 



D(l, 11) = 2 

D(lll, IV) = 2 



D(l, 11) = 2 
D(lll, IV) = 2 



D(l,ll) = 2 
D(lll, IV) = 2 



D(l, 11) = 2 
D(iii, rv) = 2 



Some Hamming distances and vectors 
of codes of length 8. 



TABLE 4.1 
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An examination of the last two examples in the table shows that 
although the (I, II) and (III, IV) pairs are not complementary their 
Hamming distance is still n/4 when compared to the original codes 
of length eight. That this feature is true in general will now be proven 
by examining and exhausting all the known ways of generating composite 
codes . 

Theore m 4.6 

Every composite code of length n (generated by one of the standard 
methods) when expressed in the sequence quadruple representation will 
have D(I,n) = D{III, IV) = n/4. 

There are presently known four general ways of generating com- 
posite codes from shorter cqdes, (there is also the method applicable 
only to codes of length 2^ given at the end of the proof). These are: 

1. Time Sequence 

51 = AB 

5 2 = AB 



2. Interlace 

T 1 ““ a 1 b -| a b . . . a ^ b 



■1"1""2"2- 



n n 



T3 = 

3. Time Sequence Exponential 
Cl Cp C dj do 

Ui = A A 



B 



m 



d d 1 di c c 1 Cl 
T T A A nn— 1 A ^ R ^ R ^ 

U2 = A a ....a B B .e«B 



Interlace Exponential 

c, d, C-, d^ c d 

Vi = A *B _‘a S -"B ™ 

d c d ,c 1 d,c, 

V2 = A b a B .....AB 



Throughout this proof it will be considered that (A,B) are a comple- 
mentary pair of length r and (C,D) are a complementary pair of 
length m. 
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<1 



= 1 is a necessary 



The parity check a. ® a , , . ® b. ® b , , . 

1 n+l-i 1 n+l-i 

condition for a sequence pair to be complementary. One of the simplest 
methods for making the parity check is to fold both codes in the middle 

9 

and double them back on themselves. This procedure is shown in 
Figure 4.1. Each 

^n/2-l^n/2 



a a 1 . . . . 
n n- 1 


n-hl - 1 


• ' • “ '^n/2+l^n/2+l 




b. 

1 


* • * ” ^n/l- l\/Z 


b b - 1 

n n 


n+l-i 


• • *°'^n/2+2'^n/2+l 



Foldover Parity Check Method 
Figure 4.'1 

column in Figure 4. 1 satisfies the parity check. If code pairs are 
written in the sequence quadruple form such as in Figure 4. 2 , it is 
seen that the columns of this array are the same as the column in the 
first array with rearrangements in order provided it is allowed that 
the columns 



a 

r 

a 

u 

b 

r 

b 



a 

u 

and a 

r 

b 

u 

b 



u r 

are the same. Since modulo 2 addition is the operation used in 
making the parity check and in determining Hamming distance, and the 
use of each column involves only modulo 2 addition which is commu- 
tative, the two columns can be considered the same. Therefore the 
fold-over method of grouping and the sequence quadruple form of grouping 
have invariant columns, order being ignored. 
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^n-3 ^n- 1 

^n^n-2 

b,b ....b ,b , 

13 n- 3 n- 1 

^n^n-2 ^4 ^2 

Sequence Quadruple Form 
of Parity Check 

FIGURE 4. 2 

Now consider the composite generating methods. 

1. Time Sequence 

Si = AB 

Sg =AS* writing in the fold- over pa rity check form . 

1 ^ J r-2 r- i r , -V /T 

b b , b, b, b, 

r r- 1 3 2 1 

a-ia^a^.e... a ^a .a 
12 3 r-2 r- 1 r f^/a sr\ - /9 

66 , 6, b, 6, D(A.B) = r/2 

r r- 1 321 

The top pair are complementary and therefore have a distance of 
r/2, similarly the bottom two rows are complementary and have a dis- 
tance of r/2 by Theorem 4. 1. To write the time sequence code in 
sequence quadruple form would not change the distance since it would be 
purely a reordering of the parity check columns^ Therefore D(I,II) 
would equal r/2 as would 0(111, IV). In terms of the composite code 
of length n D(I, II)=D(III, IV)= r/2 = n/4 whether or not (I, II) and 
(III, TV) happen to form complementary pairs, 

2. For the interlace scheme the (I, II) pair and the (III, IV) pair 
are obviously complementary sequences by construction and will have a 
distance of r/2 = n/4. 

3. Writing the time sequence exponential pair in fold-over form 

Ui=a^*a‘'^..... a'"' 

d d j where the symbol B means B time 

g m^ rev^ersed. ” 
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...a:-'- 



Us=;A "^A 


. . . A 


^1 ^2 

B B 


( 

. . . B 



1 

m 



The B's are all time reversed but the time reversal of one of the 
complementary pairs does not change the complementary property 
D(A, B)=D{A, B) = r/2. For each bit of the C code and D code there 
is a distance r/2 since every A is matched with some transformation 
of B for each of these bits. There are a total of m bits in the C 
code and in the D code, therefore the total distance of Ui folded back 
on itself is m x r/2 = mr/2. Similar reasoning gives Ug folded back 
on itself a distance of mr/2. A rearrangement of the columns will not 
change the Hamming distances Therefore D(I,Jl)=mr/2 and 
D(III, IV) = mr /2, since the tot^l length of the code is n=2mr, D(I,II) = 
D(III,IV)=n/4. 

The interlace exponential scheme is the last general method of 



code generation to be considered. Written in the fold-over form it gives 



^2 *^2 

Vi=A B A B . 



d c d , 
B >^A "^B 

d c d , 

V2=a_”^b 

^1 ^1 ^2 
BAB 



g^m/2+l^^m/2 + l 
.^m/2 + l_^m/2+l 

A ID 

g^m/2 ^^m/2 



Again each A is matched with some transform of B, D( A, B)=D( A, B)= 
D(A,B), etc. = r/2. As before, rearranging columns into sequence 
quadruple form will not effect the Hamming distances. The distance of 
the top two rows is mr/2 and the bottom two rows (Vg) also is equal 
to mr/2. The total length of the code is 2mr, therefore 
D(I, II)=P(III, rV)=mr/2=n/4. 

5. When composite codes are of length , they are sometimes 



formed in a special way from the interlaced sections of a complementary 
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pair. These sections or pieces of the codes must be of length p=2 
where j<r. Using the symbols • ^p* ^2p"^p+l^p+2‘ 

a^p and so, the following pair would be complementary if A and B 
formed a complementary pair. 



S, =A B A, B, . . . 
IP P P ^P 


np/p np/p 


=A B A^ B^ 

2p p p 2p 2p 


“^np/p^np/p 



Rather than try to use the general term to prove that D(1, 1I)=D(I1I, IV) 

= n/4 as was done in the previous cases, a short general example of 
length 8 will be used. These two codes (A,B) of length eight will be 
interlaced four bits at a time (p=2) to form a code of length 16. 

^ 24 “^ 1 ^ 2 ^ 3 ^ 4 ^ 1 ^ 2 ^ 3 ^ 4 ^ 5 ^ 6 ^ 7 ^ 8 ^ 5 ^ 6 ^ 7^8 

Writing this composite pair in fold-over form yields 

^ 14 “^ 1 ^ 2 ^ 3 ^ 4 ^ 1 ^ 2 ^ 3^4 

^8^7^6^5^8^7^6^5 

^ 24 “^ 1 ^ 2 ^ 3 ^ 4 ^ 1 ^ 2 ^ 3^4 
^8^7^6^5^8^7^6^5 

D(A,B)=4 by Theorem 4. 1. folded over is just a column 

rearrangement of the complementary pair (A,B), and since a re- 
arrangement of columns does not effect Hamming distance, the Hamming 
distance of folded-over is also 4. Similar reasoning holds for ^24. 

Since the standard sequence quadruple form is again just a re- 
ordering of the columns, the distance of the (1,11) pair and the 
(III,rV) pair is rvj A. For this example. 



I=^1^3^1^3^5^7V7 

II=^8^6^8^6^4^2^4^2 



D{I,n)=4 
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III=a j^a^bj^b^a 



D(III, IV)=4 






D(I®II,III®IV)=8. 



Therefore any composite complementary sequence pair when written in 
sequence quadruple form has D( I, II)=D(III, IV)=n /4, 



^ 

Although there is no proof that the 4 standard generating methods, 

r 

plus the special one applicable only to 2 codes, are the only way to 
form complementary codes which are composite, an exhaustive search 
of code lengths 16 and 20 revealed no codes which were not formed 
by these standard generating methods. 

The full usefulness of the* Hamming distance property is, of course, 
not yet known. The property of the distance of the complementary pair 
being equal to half the code length as shown in Theorem 4, 1 , is used 
in many of the proofs in the chapter on Hamming vectors. This property, 
D(A,B)=n/2, indicates that the code pairs are orthogonal* Because of 
this orthogonality property the leakage of an A code carrier into the 
B code filter, or vice versa, will have no effect on the output at exact 
match. 

The primary use the writer has found for the Hamming distance of 

the (I, II) and (III, IV) pairs, has been the error check that this invariant 

property offered in decomposing codes by hand. One possible use for the 

invariance of D(I,II) and D(III,IV) could be in the search for kernels of 

new lengths. If a characteristic Hamming distance for unknown kernels 
could in some way be deduced, based only on length, this invariant prop- 
erty of distance would act as a most powerful screen in the search* The 
(I, II) and (III, IV) pairs are themselves orthogonal, unless the code is a 
kernel. This may provide other possible methods for searching for new 
codes. It is difficult to say in advance what utilization will be made of any 
invariant property, but it is fairly safe to say that the property will find use. 
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CHAPTER V 

THE HAMMING VECTORS OF COMPLEMENTARY SEQUENCES 

The Hamming distance between the code pairs in a complementary 
sequence was shown by Theorem 4. 1 to be one half the length of the 
code. This Hamming distance was obtained by first finding what I shall 
call the Hamming vector, H(A,B)= A ® B or for convenience just 
H=A®B, and then counting the number of ones H contains. The Ham- 
ming vector is itself useful in recognizing the composition of a composite 
code when it is desired to break the code down into shorter lengths. This 
decomposition might be used to find previously unknown kernels, or to 
set an upper bound on the possible number of pairs of codes of any length. 

An examination of each of the methods for generating longer codes 
from shorter ones as explained in Chapter 2 will show that each method 
has a characteristic Hamming vector form, although it may be necessary 
to time reverse one of the codes to obtain this form. The method of 
generating codes of length 2^, which does not apply to complementary 
codes in general, will also be studied with its characteristic Hamming 
vector. 

In the following proofs (A,B) form a complementary code pair of 
length n while (C,D) form a complementary code pair of length m. 

The symbol (0)^ or just 0^ will mean a string of n zeros, similarly 
(01)^ will symbolize a string of n zero ones, and {a^®b^)^ will be a 
string of n elements equal to (a^®b^ . 

Theorem 5. 1 

The time sequence form of generating a composite code has the 
characteristic Hamming vector H=(0)^(1)^« 

Use the time sequence generating method on (A,B) to form 

Si= AB 

Ss= AB . ^2 



(5. 1) 



Take the Hamming vector of (Si, 82)1 H=Si ® 82 = 0^1^, this is 
seen to be true by inspection and definition of modulo 2 sum, 

jjc 

Theorem 5, 2 

The interlace form of generating a composite code has the charac- 
teristic Hamming vector H=(01)^. 

Use the interlace generating method on (A,B) to form 



Ti-a^b^a^b^a^b^. 
T2=ai bi a^b^a^b. 



a b 
n n 



a B . 
n n 



Take the Hamming vector of (Ti, T 2 ), H=Ti ® T 2 =( 01 )’^, (5. 2) 

this is seen to be true by inspection and definition of modulo 2 sum. 

>J; 

Theorem 5. 3 

The time sequence exponential form of generating a composite code 

has the characteristic Hamming vector H=(c,®d (c^®d , ... 

i m ^ m- i 

(c^®d f (d . . . d ®c, )" . 

m i i m m i 

Use the time sequence exponential method on (A,B) and (C,D) 
to form 



c 1 c^ 
Ui=A A 



^ 

.A B ^B 



B 



B 



m 



1 



d d , di c c 1 

tt a m -1 X 1 t 3 ni-n m -1 

U2=A a A B B 

Take the Hamming vector of (Ui , U 2 ), 

H=Ui® Us =(cj®d^)” •“* 

(d ®S,)". 

m 1 



(5.3) 



This is seen to be true by inspection and definition of modulo 2 sum. 
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It is to be noted that each of theTheorems 5.1,5. 2, 5. 3 have 

Hamming vectors which are anti - 6 :ymmetric about their center. This is 

required to satisfy the fold- over method of making the parity check. 
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Theorem 5. 3 also shows that along with the anti-symmetiic ^characteristic 

the Hamming vector for time sequence exponential composition has unit 

clusters of n zeros or n ones and is also anti-symmetric in both halves. 

This is easily seen because the first sum is (c , ®d ) while the last sum 

to operate on the A code is known through the parity 

check that c,®d ®c ®d, = 1 . Therefore one of the sums (c,®d ) or 

(c ®di ) must be zero and the other must be one. To summarize, for 
'ml 

the time sequence exponential form of composition the first half of the 
Hamming vector is anti-symmetric about its center (the quarter length 

point) and the first half and the last half of the vector are also anti-sym^ 
metric. 

Theorem 5,4 will concern the interlace exponential method of 
forming composite codes and it will be seen that this method will again 
form a Hamming vector anti-s ymmetric about its center, but with a dif- 
ferent arrangement of the clusters of zeros and ones. 

Theorem 5,4 

The interlace exponential method of generating a composite code has 

the characteristic Hamming vecto r H={c,®d )^{d,®c )^{c- 5 ®d , . 

^ 1 m 1 m c m- 1 ' 

(c ®d,Ad 

m 1 ' ' m 1 

Use the interlace exponential method on (A,B) and (C^D) to form 



^2 4 

Vi=A B A B l... 


c d 

, . A *^B 


d ; d , 

V^=A ""a . . 


di Cj 

, . A B . 



Take the Hamming vector of (Vi, Vg) 

H=Vi ® V 2 = (c 2®"^m- 1 1 

This is seen to be true by inspection and definition of modulo 2 sum, 

>!C 5{C 

The interlace exponential method again gives a Hamming vector 

which is anti-symmetric about its center due to the condition of satisfying 
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the fold-over form of the parity check, but its clusters of zeros or ones 

are now 2n in length rather than just n as was the case for the time 

sequence exponential method of generation. This characteristic of 

clusters of 2n in length is caused by the necessary condition of the 

parity check for the (C,D) code pair. For example, the first two clusters 

are (c.®d and (d,®c and the second two clusters are (c-.®d , 

i m i m c m- 1 

and (^2®^m 1^^* parity check to the i=l and i=2 bits 

in the (C,D) pair gives the following two equations: 

c,®c ®d,®d =1 and c^®c i®d^®d i=l . 

1 m 1 m 2 m-1 2 m- 1 

Adding one modulo 2 to both sides of both equations and rearranging 

gives (c^®d^) ®(d^®c^§l) =0 and 1 ^ 1 ® ^ 

true in general that x ® 1 = x; the refore (c , ®d ) ®(di®c )=0 and 

{c^®d 1 ) ® (d^®c , ) = 0 . The sums (c, ® d ) and (d, ® c ) must 

2 m-1 2 m-1 1 m' 1 m 

be alike to satisfy this constraint with each of the sums giving n ones 
if unlike, or n zeros if alike, and since the pairs are in time sequence 
the clusters must be in groups of 2n. A similar reasoning will show 
that the parity check for i = 2 gives identical results for the third and 
fourth terms in the Hamming vector for the interlaced exponential form 
of generating composite codes„ Therefore equation 5,4 can be written 
H(Vi. V3) = (cj ® d^)^" (""Z ® 

Before going into the special case applicable only to codes of length 

2^^ it might prove helpful to give an example of generating a code by one 

method and then decomposing it, to show that it could have been composed 

by a second method. The time sequences exponential scheme will be the 

generating method used, 

A=10 C=1001010001 

B=ll D=10000001 10 

Ui= 100101 1001 10010101 1011000000000000111100 

U2=01 1010010101010101 10001111110011001 11100 

H=Ui®U2 =1111111100110000000011 1 1 11 1 1001 100000000 
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In passing it should be noted that the shortest clusters of zeros and ones 
in H are of length 2 and that H is anti- symmetric about its center. 

If Ui is now time reversed, indicated by , the complementary 
pair is now: 

Ui=001 11 10000000000001 101 10101001 1001 101 001 
U2=01 1010010101010101 1000111111001100111100 
H = « Us = 0101010101010101010101010101 010101010101 = (01)^® 

This form of the Hamming vector indicates that H is composed of an 
interlace pair of length 20. Decomposing both and Ug by assum- 

ing an interlace pair gives 

=01 1000000101 1 10101 10 
01100000011000101001 
H = 00000000001 1111 1111 I = O^^l^® 

U2= 01 10000001011 10101 10 
10011 11 11001110101 10 
H = 111 11 11 1110000000000 = 1^°0^° 

The Hamming vector of indicates that this code of length 20 

was composed by the time sequence method from codes of length 10. 

If the form AB is assumed for the codes of length 10 then 

A=0110000001 
B=01110101 10 

This (A,B) pair is seen to be the transformation ETC^= EG on the 
original kernel of length 10 which was used to generate the code by a 
different method. Similarly the pair from Ug is also a transformation 
of the original (C , D) pair« Therefore the given code of length 40 
could have been generated by two different methods. 

A second example will now be given using the same (A,B) and 
(C,D) kernels as were used in the first example, but his time utilizing 
the exponential interlace scheme of generation rather than the time se- 
quence exponential form. 
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A=10 

B=ll 



C=1001010001 
D=10000001 10 



Vi=10 1101 0001 00100001 001 00001 0001 1101 111000 
Vs=01001011 10 1101 11010001 11010001 1101111000 
H=Vi®Vs =111111111111111100001 1110000000000000000 

Note that the smallest cluster is four, as would be expected of an inter- 
lace exponential scheme with the (A,B) pair of length two. The H 
vector is again anti- symmetric about its center. Time reversing Vi 
give s 

V^=000 11 11011 10001 00001001 00001001000101101 

v7=01 001 01 1101 1011 10100011 10100011 101 111000 

20 

H=Vj^®V2 =0101010101010101010101010101010101010101 = (01) 

As before this form of H indicates an interlace code pair. Decomposing 

both V ^ and into interlace code pairs gives 

=001 11101000100010110 
~ 01101000010001000011 
H =01010101010101010101 = (01)^° 

V2= 001 11101000100010110 
10110111101110111100 
H = 10101010101010101010= (10)^° 

Vj and Va both break down into interlaced pairs of length 10„ The top 
row of breaks into A=01 10000001 , B = 01 110101 10 which is the 

transfornnation EG operating on the original (C^D) pair as in the 
previous exannple. Each of the codes for Vs would break down in 
similar fashion to some transformation of the original (C,D) code pair. 

One use for the method just shown could be to search for larger 
kernels than those currently known. An example might be kernels of 
length fifty found through the use of time sequence exponential or inter- 
lace exponential, where both (Aj,B) and (C,D) are codes of length 10. 
These would form codes of length 200^ From this as a starting point, 

a decomposition to codes of length 100 and then to length 50 might 
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I. 



be possible through this technique^ This technique will be discussed 
much more thoroughly towards the end of this chapter, but it is neces- 
sary to consider the special case of 2^ code composition and the 
Hamming vectors of the sequence quadruple form before continuing with 
this interesting possibility. 

The special generating method which is applicable only to codes of 

length 2 is actually a generalization of both the time sequence and the 

interlace methods. As mentioned in the previous chapter, this method 

takes a code pair of length 2 and from them forms a code pair of 
r Hh 1 

length 2 , The allowable combinations are formed by interlacing 

pieces of the A code and the B code of length 2^, where m=, 0, 1 , 2, . ,r. 
As an example a code of ^length 8 can be used to generate codes of 
length 16, In this case r=3 and m has the possible values 0,1,2, 3, 
These values of m give section lengths of 1,2, 4, 8 respectively, 

A— a a 2 a ^ a^ a ^a ^ a ^a g 

for m=0, 

S21=^1^1^2^2'"3h^4V5^5^6^6^7^7^8^8 

=0101010101010101 =(01)® = (oSV 

This is exactly the same as the interlace procedure which is a standard 
generating method. 

The section lengths are 2 bits long for m=l, 

^12“^1^2^1^2^3^4^3^4^5^6^5^6^7^8^7^8 
^ 2 2 1 ^ 2 ^ 1 ^ 2® 3 ^4 ^ 3 ^ 4 ^ 5^ 6^ 5^ 6^ 7^ 8^ 7^ 8 

H=S^2®S22 =00110011001*0011 = (O^l^)"^, 
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m=Z gives 



^14 

S24=^1^2^3^4^1^2^3V5^6^7^8^5^6^7^8 

H=Sj^®S24 =0000111100001111 = (0^1^)^, 

For m=3 the lengths are eight, which is the entire code. 

S,g=AB 

^28=^®- 

This is the same as the time sequence form giving a Hamming vector 
H=0000000011111111 = (0®1®). 



The Hamming vectors from this special construction are altering 
clusters of zeros and ones, each cluster being the length of the segment 
used in the generation. 

Theorem 5. 5 

r-f* 1 

Complementary sequences of length 2 formed from an (A,B) 
r 

pair of length n=2 using the special generating method with segments 

2n/ m 

of length m have a characteristic Hamming vector H=(0 1 ) 

Use the special generating method utilizing segments of length m 
of the pair ( A, B) to form 

^lm"^1^2^3” * ‘^m ^l’^2* ‘ *^m ^m+1®' m+2. . .^2m ^m + l^m+2* ‘ 
^2m~^1^2^3* • ‘^m ^1^2’ * “^m ^m+l^m+2* ’ ’ ^2m ^m+l^m+2* * * * 



Take the Hamming vector of ^2m^ 

H=S, dlS, = 

Im 2m 

This is seen to be true by inspection and definition of modulo 2 addition. 



ij< Jlc >Jc 

It now seems suitable for completeness to include a theorem which 
has been alluded to many times but not proven. This is the anti- sym- 
metric property of all Hamming vectors of complementary sequences. 
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Theorem 5. 6 

The Hamming vector of a complementary sequence pair has a pat- 



tern which is anti- symmetric about the center of the Hamming vector. 



1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 



•A." a. a-ja-^a. 
12 3 4 



a , a 
n- 1 n 



B=b,b^b-b . 
12 3 4 



b ,b 
n- 1 n 



H(A,B)= h^h^h^h^. . . h^^2* • « \ 



a. ® a . , .§ b. ® b . = 1 

1 n+ 1-1 1 n+l-i 



is a complementary pair of 
length n 

by definition of Hamming vector 
by equation 2. 4 



a. ® b. = h. 

1 1 1 and 

^n+l-i ® ^n+I-i “ ^n+l-i 



both by definition of h. 

^ 1 



Add ^ ® 1 ) to both sides of the equation 

h , . = h. 
n+ 1-11 * 



Make a change of reference to the center of the vector, which lies 
between m/2 and m/2 + 1, rather than the end by having i=n/2 -r then 
n + l-i=n/2+l + r 



9. h = h +l+r for all 0 < r <n/2 -1 by step 8. 

n/2-rn/2 —— / jr 

1 0. Therefore the Hamming vector is anti- symmetric about its 
center when it is formed from a complementary pair. 



The Hamming vectors of (I, II) and (III, IV) in the sequence quad- 
ruple form are invariant under transformation if H(I,II) and H(III, IV) 
are taken as an unordered pair. The following four theorems which 
parallel Theorems 4. 2 , 4. 3 , 4. 4 , 4. 5 on Hamming distances will dem- 
onstrate this characteristic of unordered Hamming vector pairs. 

Theorem 5. 7 

A change in the order of (I, II) or (III, IV) or both will not change 



H(I,II) or H(III,IV). 



70 



1. H(U, V)=u^ ® by definition. 

2. u. ® V. = V. ® u. since modulo 2 addition is commutative, 

1111 

3. therefore H{U,V) = H(V,U) 

>•< 



Theorem 5. 8 

Complementing the pair (I, II) or (III, IV) or both will not change 
H{I,II) or H(III,IV). 



1. 


H{U, V) = u. « V. 

11 


by definition. 


2. 


U. ® V. = U. ® V. 
1111 


by Theorem 2. 1 , 


3, 


therefore H{U,V) = 


= H(U,V). 
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In the following proof it is assumed that each vector I, II, III, IV is 
of length m. It is also necessary to show that x = x. This is easily 
seen sincex® l=x and x ® l = x, therefore x = x ® (1 ® 1) = x. 



Theorem 5. 9 

Complementing one of the pair (I, II) and one of the pair (III, IV) 

exchanges the Hamming vectors H(I,II) and H{III, IV). 

I— a-ia^ar-.......a -• 

13 5 n- 1 

II=a a a^ 

n n- 2 2 

III=b, b b,. ...... b . 

L 5 D m- 1 

IV=b b -2 b^ 

mm 2 

1. Hd.ID^hjh^hjh^ 

2. Hdll.IVj^hjh^hjh^ 5i^ 

cause of the parity check. 

3. H{1, II)=H(I,IT) =hjh2h3h^. . h^. 

4. H(m, rv) rHClII,"!?) = =hjh^h 2 h^ h^ , 

5. Therefore H(I, II)=H(m, IV) =H{IU,IV), 

6. and H(III, IV) = H(I, II) =H(1, 11], 

Jie >:< « 
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Theorem 5, 10 



Exchanging the (I, II) pair with the {IU,IV) pair exchanges their 
Hamming vectors. 

1. A change in the Hamming vector can only be caused by a bit 
change. 

2. Exchanging the (I, II) pair with the (III, IV) pair changes 
position but not bits* 



Earlier in this chapter codes of length 40 generated from the quad 
were decomposed into kernels of length 10 by use of the Hamming 
vector. The conjecture was made at that point that perhaps codes of 
length 200 manufactured from kernels of length 10 could be decomposed 
into kernels of length 50. In order to prove this decomposition is not 
possible several more theorems are necessary. These theorems will 
be proved before preceding to the decomposition problem. 

Theorem 5.11 

Under all possible general transformations of a complementary pair 
where n > 2 there are only four possible Hamming vectors, two of 
which are complements of the other two. 

1. (A,B) are a complementary pair of length n where n > 2, 

2. H(A, B) = H interlaced with H , where H are the odd bits 

of the Hamming vector and H^ are the even bits of the Hamming vector 

in reverse order. The symbol H=H^ H^ will be used for H^ inter- 
laced with H in exactly the same manner that (I, II) and (III, IV) are 
e 

interlaced to form A,B. Also h^ = a^®b^ as before. 

3. Let Hj = (I»III)=a^®bj^, a^»b 3 . .. = h^h^h^.. . h^_ ^ , 

H2=(II»IV)=a^»b^.a^_2»b^_2...a2«^b2 = Vn-2'-^2' 

H2=(I® IV)=a^®b^,aj©b^_2- .... ^n- ’ ’ ^n- 1 ’ 
H4=(II®III)=a^«b^ . a^_ 2®b3. .... a2®b^_ j 
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by Theorem 2. 1 , 



4, u.® V = u.® V. 

11 11 

therefore 

(r®TTi)=H^ (ir«Tv)=H2 (i<i)iv)=H^ (mm)=H^. 

5* If h. = u.®v., then h.® 1= h. = u.®v. ® 1 = u.® v. = u.® v. 

Ill 1 111 1111 

therefore 

H(U,V) = H(U, V)=(U,V) 

and H^=(T®III) = (I®iTj) H^=(T®IV) = (I®lV) 

H2=(Ti»iv) = (II®IV) H^=(Tr®iii) = (Ii®m). 

6. = h^h^h^ h^. 

* H, = h h , h,, 

2 1 n n- 1 1 

7. But h /9 =b /^ , 1 , by Theorem 5.6. 

n/2- r n/Z + 1 + r ^ 

8. ^n/2" ' • ^2^ 1 

H 2 * Hi = hi ••”^n/2 ^n/2'**^2^1’ 

9. Therefore Hi * - HJ^“H^ = H^ * H i* 

10. H, >!« H. = K.K. ..... t . 

3 4 i Z n 

11. H^ H^ = 1 .... . 1 ^ 1 . 

^n/2- r"^n/2 + l+r. 

13. Therefore H^ * H^ = = H^ * H^. 

14. Table 5. 1 is an exhaustive list of transformations on a 
complementary pair in the sequence quadruple form with the Hamming 
vector for each transformation. 

and their complements are the 

only possibilities which exist. Therefore only four Hamming vectors 
are possible and by steps 9 and 13 they occur in pairs which are 
complements of one another. 

>!< ^ ❖ 
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TRANSFORMATION 



I 


(I II 


III 


IV) 


Cl 


(T 


ii 


III 


IV) 


Cz 


(I 


II 


Hi 


W) 


c 


(T 


Ti 


TTi 


TV) 


^1 


a 


II 


Tii 


IV) 


^2 


(I 


IT 


III 


iv) 


^1 


{II 


I 


III 


IV) 


Tz 


(I 


II 


IV 


III) 


T 


(II 


I 


IV 


III) 


Z 


(II 


1 


TiT 


'IV) 


Y 


(IT 


■ I 


TTi 


IV) 


X 


(H 


■ I 


III 


W) 


w 


(II 


1 


III 


Tv) 


V 


(T 


II 


IV 


TTi) 


u 


(T 


II 


Tv 


III) 


s 


(I 


TT 


■ Tv 


III) 


R 


(I 


IT 


IV 


III) 


Q 


(II 


T 


IV 


in) 


P 


(n 


■ I 


Tv 


III) 


O 


(Ti 


1 


III 


IV) 


N 


(II 


I 


TTi 


TV) 


M 


(T 


ii 


IV 


III) 


L 


(Ti 


1 


Tii 


iv) 


K 


(T 


Ti 


Tv 


TTi) 


J 


(n" 


■ T 


Tv 


Hi) 


H 


(I 


II 


Tv 


TTi) 


G 


(n 


■ T 


IV 


III) 


F 


(I 


Ti 


TTi 


IV) 


D 


(T 


II 


III 


fv) 


tr 


(II 


I 


Tv 


TiiF 


B 


(H 


■ I 


IV 


Hi) 


0 


(II 


1 


Tv 


III) 



HAMMING VECTOR 
H{I)=Hj’^H^ 

H(C) 

H(A^)=H^*H^ 

H{Tj)=H^*H3 

H(T2)=H3*H^ 

H(T) =H3=i'Hj 
H(Z) =H4*H3 
H(Y) =H4*!‘H3 
H(X) =H^*H3 
H(W) =H4’^H3 
H(V) =H3*H4 
. H(U) =H3*H^ 

H(S) =H3*H4 
H(R) =H3*H4 
H(Q) 

H(P) 

H(0) =H4*H3 
H{N) =H4*H3 
H(M)= H3*H^ 

H{L) =H^*H3 
H(K) =H3«H4 
H(J) 

H(H) =H3*H4 
H(G) 

H(F) =H^*H3 
H(D) =H^>1'H3 
H(ir) =H 2 *H^ 
H(B)=H3*H^ 
H(0)=H^*Hj 



COMMON VECTOR 
H(I) 

H(I) 

H(I) 

H(I) 

H(I) 

H{I) 

H(T^) 

H(T^) 

H(I) 

H(T^) 

H(T^) 

H(T^) 

H(Tj) 

H(T^) 

H(T^) 

H(T^) 

H(T^) 

H(I) 

H(I) 

H(T^) 

H(T^) 

H(T^) 

H(T^) 

H(T^) 

H(I) 

H(T^) 

H(I) 

H(I) 

H(I) 

H(I) 

H{I) 

H(I) 



The E transformation just permutes the first two and last two pairs 
of (I II III rV) and since modulo 2 addition is commutative no change 
would take place in the H vector. The last column identifies all 
Hamming vectors in terms of two vectors and their complements. 

TABLE 5. 1 
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Theorem 5.12 



If a complementary pair ( A, B) of length n is operated upon by a 
complementary pair (C , D) of length m to form codes Ui, Ug by the 
time sequence exponential method, or to form codes , Vg by the inter- 
lace exponential method and either Ui or U 2 is time reversed, symbo- 
lized by U or if Vi or V 2 is time reversed symbolized by V then 

H{U) = H(V) = H(A^ , ^ b^) pair indicates some 

general transformation of the original (A,B) pair. 

1. A ® B = H(M) . where B indicates the time reverse of 

“ the original B 

2. A « B = H(H). 

3. B«A = A®B =H(Ti). 

4. B ® A = A ® B^H(L). 

5. Referring to Table 5. 1, 



H(M)=H(Ti), 
H(H)=H(Ti), 
H(B) = H(Ti), 

^1 ^2 

Ui=A A .. . 
d d 

Ua=A 



c d, 
A B_ *. 

d 1 c 
A ' B "! 



B m 



B 



and 
Vi=A 



/i 



c d 
m B m 



d c d , d, r, 
Vs=A "" B ” A “:.‘.A * B 



Time reversing one of the pair in each, for example, the use 
of T 2 yields 



Ui = A A . . 



Us = B 



1 



B 



and 



1 ^2 
Vi = A B ^ A 



c 

c 



1 "^1 '“2 
V2 = B A B . 



c d, 
A_"^B \ 
c d, 
B ”^A ^ 



B 



m 



m 



c d 
A ”^B 



B "^A 
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7. 



H(U) and H(V) consist of clusters of (A ^ ® B and 
“ d. d. 

clusters of (B ^ ® A ^). 

8, The Hamming vectors for these clusters have the following 
four possible forms: 

A® B,A®B, B®A, and B ® A. 

These four possible vectors were seen in steps 1,2, 3,4, 5 to all be 
equal to H(Ti). If Ui and Vi had been time reversed, using the 
operation Ti, the result would have been H (Ti), 

1 0. Since all the clusters are identical and there is a cluster for 
each of the m bits in C and each of the m bits in D, 

H(U) = HV) = [H(Ti)]^,^ for Ug or Vg time reversed or 
in general 

H(U) = H(V) = H{A^,B^) where (A ,B ) signifies some transfor- 

mation of (A ,B). 

>;c >;< 5*^ 

Theorem 13 

The kernel of length 2 when operated on by a complementary pair 
by either the time sequence exponential or interlace exponential method 
can always be decomposed at least once by the interlace method. 

1. If one of the code pair formed by either exponential method 

is time reversed, H(l^) = H(V) = H(ASB^) by Theorem 5. 14. 

2. There are only two possible Hamming vectors for the 
quad, 01 and 10. 

3. Therefore H(U-) = H(V< ) = (01)^^ or (10)^^ both of which 
are decomposable by the interlace method. 

>!c >;< 
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Theorem 5, 14 



A complementary pair (Ui, Ug) formed by the time sequence 
exponential method cannot be decomposed by the time sequence method. 

1. Let (Ui , U 2 ) be a complementary pair of length 2mn formed 
from the complementary pair (A, B) of length n operated 
upon exponentially by the complementary pair (C,D) of 
length m, formed in the same manner as in Theorem 5. 3. 

2 . H(U)={c, ® dm)’^(c.»d i)"...(c ®dj"{d,»a )"...{d ®c , 

3. The first mn bits in H(U) are m clusters composed of n 
zeros or n ones dependent upon the bits of H{C,D). 

4. By Theorem 5. 1 H=0 1 or 1 0 , for decomposition 

of a code of length. 2mn by the time sequence method, 

5. This requires the distance of (C,D) be zero or m, but by 
Theorem 4.1 D(C,D) = m/2, therefore H(U) is not decom- 
posable by the time sequence method. 

6. If Ui or Ug is time reversed H(U) = H(A^,B^)^^ by 
Theorem 5.12. 

7. This is incompatible with the form since it would 

require the first m/2 clusters of H(A^,B^) to have a distance 
of zero and the second m/2 clusters of H(A^,B^) to have a 
distance of n. Both are impossible by Theorem 4, 1 which 
requires a distance of n/2. 

>!« ❖ 

Theorem 5.15 

A complementary code pair (VijVs) formed by the interlace ex- 
potential method (using the same definitions and symbols as Theorem 5. 4) 

can be decomposed by the time sequence method only if H(C,D)= 
QmyZ^iWZ^^ ^m/2 Qtn/Z. 
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1. H(V) = (c «d )^"(c,«d (c €d,)^" by Theorem 5. 4 

1 m Z m- 1 ml 

and equation 5. 5, 

Z. The clusters of zeros and ones in this Hamming vector are in 

lengths of Zn and are dependent in value upon the bits of H(C,D), 

3, is the requirement for decomposition of this code 
by the time sequence method. 

4, If H(C , = or its complement, it will satisfy step 3 

since (0 ) (1 ) =0 1 =H(V), and it also satisfies 

Theorem 4.1. 

5, If Vi or Vs is time reversed H(V) =H(A^ , by 

Theorem 5. I Z. 

6, This is incompatible with the form H=0^^l^^ since it would 

require the first m/Z clusters of H(A^, ) to have a distance 

of zero, and the second m/Z clusters to have a distance of n. 
Both are impossible by Theorem 4. 1 since it requires a 
distance of n/Z. 

7, Step 4 showed that H(C , D) = 0^^^1^^. This satisfies the re- 
quirement for the time sequence decomposition by inspection. 

Any other arrangement of H(Cj,D) would fail to allow decom- 
position. 

>JC 



Theorem 5.16 

Any complementary sequence pair (Vi, Vs) formed by the interlace 
exponential method as in Theorem 5.4 (using the same definitions and 
symbols), the quad as the (C,D) pair, is time sequence decomposable. 

I, The quad has only two possible Hamming vectors, 01 and 10. 

Z. This satisfies the form H=0^^^l^^^ where m=Z, the quad 
length. 

3. Vi and Vg are therefore time sequence decomposable by 
Theorem 5.15. 
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Theorem 5.17 



A complementary pair formed by either exponential method is de- 
n/2 

composable by the interlace method if and only if H(A^, B^) = (01 ) or 

n / 2 

(10) . All definitions and symbols are assumed to be the same as 

those in Theorems 5. 3, 5. 4 and 5. 12. 

1. By Theorem 5. 2 these codes are decomposable by the inter- 
lace method only if they are of the form H=(01)^^ or (10)^^. 

2. H(U) and H(V) are clusters of at least n ones or n zeros 
and therefore cannot have the form (01)^^. 

3. H(U) and H(V) = H( A^ by Theorem 5.12. 

4. If and only if has the form (01)^^^ or (10)^^^ can 

H(U) or H(V)= [(01)'^^^] = (01)^'^, or using 10 instead 

of 01, (10)^^, This satisfies the requirement of Theorem 

r r 

5.2 that the Hamming vector be of the form (01) or (10) . 

>|< 

A brief description of the method used to prove that codes of length 
200 formed from kernels of length 10, cannot be decomposed into ker- 
nels of length 50, will now be given as a guide to the actual proof of the 
theorem. The method used is to exhaust all possible methods of de- 
composition by examining all possible Hamming vectors for the (A,B) 
and (C,D) pairs. Table 5.2 lists all possible Hamming vectors for 
kernels of length 10. 

Kernel 1 Kernel 2 

A= 10000001 10 A=1001 1011 11 

B= 1001 010001 B=1 100001010 

The 4 possible Hamming vectors. The 4 possible Hamming vectors. 

0001010111 0101100101 

1110101000 1010011010 

0000101111 1100101100 

1111010000 001 1010011 

All possible Hamming vectors of kernels of length 10. 

TABLE 5.2 



The two methods of generating codes of length 200 from kernels 
of length 10 are the time sequence exponential method and the interlace 
exponential method given respectively by 



^2 ^10 ^1 
Ui=A A A B . .. 


. . .B^IO 




U3=A '^A . a "b " \ . 




(5.6) 


v,=a"‘ b”*' b^... 


...a'*“b‘*‘'> 

...a">V‘ 


(5. 7) 



where (A,B) and (C,D) are kernels of length 10. 

The Hamming vectors of these codes of length 200 are anti- sym- 
metric about their centers and are in clusters of 10 zeros or 10 ones 
for H(U), and in clusters of 20 zeros or 20 ones for the H(V) vector. 
Time reversing either Ui ot Ug and similarly either Vi or Vg , 
gives 20 clusters of Hamming vectors of the forms listed in Ikble 5.2, 
since each of the 20 A codes is now matched with some transform of B. 

The known possible ways for kernels of length 50 to generate com- 
posite codes of length 20 are: 

1. Either interlace or time sequence kernels of length 50 to 
form composite codes of length 100 and then either inter- 
lace or time sequence these n=100 codes to form composite 
codes of length 200. 

2. Take kernels of length 50 and either time sequence expo- 
nential or interlace exponential with the quad to form codes 
of length 200. 

3. Take kernels of length 2 and either time sequence exponen- 
tial or interlace exponential with kernels of length 50. 

An examination of all possible Hamming vectors of length 200 from 

each of these methods will show the incompatability of these Hamming 

vectors with those generated by the kernels of length 10. 

Theorem 5.18 

No codes of length 200 formed from kernels of length 10 can be 



decomposed into kernels of length 50 by standard methods of decomposi- 
tion. 80 



Given Ui Ug and Vi Vg both complementary pairs of length 
formed by the time sequence exponential and interlace expo- 
nential methods as given in equations 5. 6 and 5« 7 respectively. 
The time sequence method of decomposition is impossible for 
each of the following Hamming vectors for the following reasons. 

a. H(U) can never be decomposed by time sequence. 

Theorem 5. 14. 

5 5 

b. H(V) requires a form 0 1 by Theorem 5. 15. 

Table 5. 2 contains no vector of this form. 

c. H(V ) and H(U) can never be decomposed by time 
sequence. Theorems 5.14,15. 

The interlace sequence method of decomposition is impossible 

for each of the vector forms for the following reasons. 

a. H(U) and H(V) can never be decomposed by interlace. 
Theo rem 5.1 7. 

b. H(V) and H(U ) are decomposable by interlace only if 
H(A^, B^) is of the form (01)^ or (10)^ byTheorem 5.17. 
Table 5. 2 contains no such vectors. 

Assume that one of the Hamming vectors H(U) or H(V) has a 

second exponential form, first with the kernel of length 2 

operated on exponentially by the kernel of length 50, and 

secondly by the kernel of length 50 operated on exponentially 

by the quad. L.et H^ be the designator for the Hamming 

vectors of the 2 and 50 combinations. 

a. Consider first the quad acted on exponentially by the 

kernel of length 50. H^U)=H^(V)= H^(A^,B^). But since 

A^,B^ is the quad, by Theorem 5. 13, it is decomposable 
by the interlace scheme. This is not possible since 
step 3 above exhausted all possible interlace decompo- 
sitions. 

b. Next consider the kernel of length 50 acted on exponen- 
tially by the quad. Considering first the exponential 
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interlace scheme, by Theorem 5. 15 if H^(C,D)=01 or 10, 
which it does for the quad. Therefore the code is time se- 
quence decomposable. However, in step 2 all of these time 
sequence possibilities were exhausted and therefore this de- 
composition is not possible. H^(U) would have a form con- 
taining either 50 zeros or 50 ones in clusters. For the 

original 10 length codes to form these would require a 
5 5 

0 1 Hamming vector which is not listed in Table 5. 2. 

This exhausts all possible Hamming vectors that might lead to de- 
corpposition, therefore codes of length 200 formed from kernels of 
length 10 cannot be decomposed into kernels of length 50 by standard 
methods. 

❖ 5jc 

Table 5.3 shows the pcfssible Hamming vectors for codes of 
length 26. 

A=0100 1101 111010111100111010 

B= 10 11 001 00001 nil 1100111010 

The 4 possible Hamming vectors. 

11111111111101000000000000 

000000000000101 11111111111 

11101110111010100010001000 

0001 000100010101 1101 1101 11 

All possible Hamming vectors of the kernel 
of length 26. 

Table 5. 3 

An examination of Table 5, 3 which is exhaustive for the Hamming 
vectors of length 26 and Table 5.2 for 10, shows that the next two 
theorems can be proved by identical methods to those used in Theorem 
5. 18. Since these proofs exactly parallel that of Theorem 5. 18 they 
will not be given. 

Theorem 5.19 

Complementary pairs of length 520 formed from the time sequence 
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exponential or interlace exponential from kernels of length 26 and 10 
cannot be decomposed by standard methods into kernels of length 130. 

5{c 5|< 



Theorem 5. 20 

Complementary pairs of length 1352 formed from the time sequence 
exponential or interlace exponential method from kernels of length 26 
cannot be decomposed by standard methods into kernels of length 338. 



The study of Hamming vectors thus far has disclosed that a pattern 
n/2 n/2 n/2 

of (01) or 0 1 or their complements in either the (A,B) 

pair or (C,D) pair, is required in order that a code formed from one 
kernel length by an exponential method be decomposed into a pair from 
a different kernel length. The quad is only kernel which thus far satis- 
fies this criteria. However, the other kernels come close, as a check 
4 4 2 2 

shows that 0 101 and (01) 10(01) are among the Hamming vectors of 

12 12 

length 10, and 0 101 is a Hamming vector of length 2 6. 

If the quad is involved in the formation of composite codes it makes 
their decomposition more likely^ This is shown in Theorem 5. 12 and 
Theorem 5. 16. Appendix II demonstrates this fact quite clearly since 
it contains all codes of length 16, less some operational redundancies. 
Each of the code pairs and its time inverse is decomposable by some 
method into shorter length codes. Appendix III lists all the pairs of 
length 20, less some operational redundancies. Each code pair or its 
time inverse can be broken into its originating pair through the use of 
Hamming vectors, but no code and its time inverse can both be decom- 
posed by the standard Hamming vector methods. 



The writer would like to end this chapter with a conjecture that 
kernels of length 50 do not exist. This conjecture is based upon three 
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bits of evidence, the proof of Theorem 5„ 18 being the first*, The second 

2 2 2 2 

piece of evidence is based upon the fact that 50=7 +1 =5 +5 , From this 
one might assume that codes of 200 formed from 10”s might cover just 
one set of the ones determined by equation 2« 5„ This proved not to be 
true however, since the number of ones in both the exponential formations 
of codes of length 200 from kernels of length 10, if they were decompo- 
sable into kernels of length 50, would exhaust all possible unordered 
pairs of ones for n=50 as determined by equation 2*, 5« The last bit of 
evidence to base this conjecture upon is that the code length is not twice 
a prime number. The only previous kernel which might have existed, but 
was not twice a prime number, was n=18, and it did not exist. Similar 
reasoning extends this conjec.ture to kernels of length 130 and 338. 
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CHAPTER VI 



SUPPLEMENTARY AND CYCLIC COMPLEMENTARY CODES 

This chapter is concerned with two classes of codes. The first of 

these classes consists of quadruples of sequences with the property that 

the total number of likes at each spacing equals the total number of un- 

2 

likes at the same spacing. These are called supplementary codes. The 
second class of codes, which will be defined later, consists of the cyclic 
complementary codes. Both of these code types have complementary 
sequence pairs as a possible subset; this allows the use of their prop- 
erties as screens in the search for new kernels. 

Supplementary codes will be discussed first, since the n=26 search 
involved the use of this property as a screen, while the cyclic comple- 
mentary property was not used until the search for kernels of length 34. 
Only one theorem will be proved for each of these types of codes, since 
these two theorems were the only ones applied in the actual search for 
new kernels. 

If u , u^u^ u 

1 Z 3 n 

V, v^ v^ V 

12 3 n 

w , w^w^ w 

1 Z 3 n 

x, x^x^ X are sequences of zeros and ones and 

12 3 n 

satisfy the constraint of equation 6. 1, they form a quadruple which is 
supplementary. 

i=j 

S (u.®u , . .®1) + (v.®v , . .®1) + w.®w . .®1) + x.®x , . .®1) = 

' 1 n+i-j ' ' 1 n-l-i-j ' 1 n+i-j i n+i-j 

i=l 

j 

Z (u.®u .) + v.®v , , .) + (w.®w . .) + x.®x . .) = 2j (6. 1) 

1 n+i-j' 1 n+l-j ' 1 n+i-j' i n+i-j' \ / 

for all j, l^j_<n - 1. Note that there are 2j like pairs and 2j unlike 

pairs at a spacing of n-j. 
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For simplicity in the proof, we use the symbolism that was used in 






Chapter 2, 
let 



Ly, = S u. ® u . . .01 
U 1 n+i-j 



J 

S 

i=l 

j 

S 

i=l 



for all j, 1 j < n- 1 



Uf, = S u. ® u , . . 
U 1 n+i-j 



A restatement of equation 6. 1 in this symbolism is 



Theorem 6. 1 

Any complementary sequence pair, (A,B), written in standard 
(1 II III IV) form has the property that I, II, III and IV are a 
supplementary quadruple. 

1. L^=Ug and by the definition of complementary. 

2. Considering just the even values of j rather than all values, 



3. 

4. 

5. 



La~Li + Ljj, 

Ua=Ui+Uii, 

^B"^III ^IV’ 

Ub=Uiii + ^IV 

Li + Lji =Uiii + Ujy. 

^iii''’^iv"^i ^ir 

Adding the equations in steps 3 and 4 gives 



Li + Lji + Lju + Lj^ - Uj + Uji + Ujii + 



6. This is the same as equation 6, 2 and I, II, III, IV formed 
in the sequence quadruple form from a complementary pair 
are supplementary. 

Since the supplementary property applies only to the even spacings 
in a code pair, the converse, that all interlaced supplementary quad- 
ruples are complementary, is not true. 

^ 
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The cyclic complementary property is useful in itself for a communi- 
cations or telemetering system; however, at this time our concern is 
for its property as a necessary condition for a code pair to be complemen- 
tary. A cyclic sequence or code, as the name suggests, is a never 
ending sequence of zeros and ones which has a period of n bits. A cyclic 
complementary sequence pair, is a pair of cyclic codes, each of period 
n, where the number of likes of one sequence equals the number of un- 
likes of the other sequence for n possible matches for all spacings 

from 1 to n-1, or stated in terms of the bits of the A and B codes 
n i=n 



C . = E a.® a 
J i=i ^ 



= E b.i 



n-j+i i n-j+i 



1 < j < n- 1 , 



(6.3) 



where a , . = a. and b = b. since the period is n. 

Theorem 6. 2 

A complementary sequence pair (A, B) of length n, if written in a 

cyclic fashion, is always a cyclic complementary pair of period n. 

1. Since (A, B) are a complementary pair they satisfy 
equation 2.3. 



3. 



4. 

5. 



J 

f. — E a.® a ... — ut i^.w ‘ 

J 1 n-j+i 1 n-j+i 



i=J 

= E b.® b 
i=l 



1 for j, l = j=n-l. 



Expanding for fj^ gives 

a,®a =b,®b®l which is equal to a®a,=b®bT®l. 
Inin nlnl 

Expanding for f^ ^ gives 

(a^® a^) + (^2® ^3^ ^ ^3® ^4^ +. . . . si^)={h^®b2®l) + 

(b^® b^® 1 ) + (b^® b^® 1 )+..... . +(b^^ 1® ^ 

Adding the equations of step 2 to step 3 gives 
Similarly expanding f^ gives 



(a, ® a 1 ) + (a^® a ) = (b , ® b^ , ® 1 ) + (b^® b ® 1 ) 
1 n-1 2 n 1 n-i n 

which is equal to 



(a^_ J® aj) + (a^® a^) = b^® 1) + (b^® b^ ® 1). 
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t: 

ir 



6 . 

7. 

8 . 



Expanding ^ gives 

(a^«a3) + (a2®a^) + (a3®a3) + , . . ^ 1) + 

(b2®b^$l) + (b3®b3®l) + + 

Adding the equations of step 5 to step 6 gives C^. 



Continuing this same procedure through all possible values of 
j would show that = f ^ + f^ y Therefore the cyclic com- 
plementary constraint equations are based on the sum of two 
restricted portions of the complementary constraint equations. 
Therefore all complementary pairs have to be cyclically 
complementary. The converse is not nece ssar y true since 
Cj is a sum, and in a sum the addend and augend are not 
unique. 



Since Cj = f ^ + f^ y each constraint of the complementary property 
is used twice as j varies from 1 to n-1; therefore Cj is symmetric 
and is centered at application of the cyclic complementary 

property in Chapter 8, half of the characteristic cyclic number is de- 
leted due to this symmetry property. 

Both of the screens described in this chapter were suggested by 
14 19 

Dr. Golay; ’ however, to the best of my knowledge the theorems 
and their proofs have not appeared in print. 
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AN EXHAUSTIVE SEARCH FOR KERNELS OF LENGTH 26 

CHAPTER VII 



In his paper Golay emphasized the need for an exhaustive search for 
complementary code pairs of length 26. This chapter describes such a 
search. The purpose of this search was to determine if there were any 
kernels of length 26 and if any of these exist, the number of such 
kernels. It would prove very difficult to surmise the number of kernels 
in a code of this length predicated on the known kernels of other lengths, 
since n=2 had one kernel, n=lO had two, and n=18 had zero. An ex- 
haustive search for kernels of length 26 might throw some light on a 
possible general method for finding kernels of longer lengths, or might 
disclose some sort of pattern showing the distribution of kernels among 
the possible code lengths. 

The only feasible method to accomplish this search was with a high 
speed digital computer. This chapter is concerned with the computer 
program and the results of the exhaustive search for kernels of length 26. 

At first glance it would appear that there are 2^^ possibilities to 
be screened for an exhaustive search for all possible kernels of length 
26. Although this is true, the application of some of the theorems de- 
veloped in the earlier chapters immediately eliminates from considera- 

52 

tion large blocks of the 2 possibilities which are either redundant or 
impossible. This is a very necessary procedure since even with the 
highest speed computers of today it would be impossible to investigate 
this nuirnber of possibilities in a life time. 

The program was developed for the GDC 1604; this computer ac- 
complishes approximately 200,000 operations per second. This speed 
coupled with the powerful screens used in programming made possible a 
reduction of the computer run-time to approximately 75 hours. A brief 
outline of the various screens used and their reduction factor will now be 
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given before continuing in more detail with the actual programming 
techniques used. The word code or code pair as used in this chapter 
and also in the next chapter will be understood as possible complemen- 
tary code pair. 

The first screen utilized was the number of ones which must appear 

2 2 

in each of the code pairs. This was given by n=(n-p-q) + (p-q) 

as derived in chapter 2, where n is the length of the code, p is the 

weight of the A code and q is the weight of the B code. The possible 

solutions for n=26 were the unordered weight pairs (16,15), (16,11), 

(11, 10), (15, 10). The pair (p=l6, q=15) was arbitrarily chosen to be 

used in the program. This reduced the total number of possible code 

pairs from 2 ;::i8xl0 down to — x ^ j ^ | ^ 4x10 . This 

was a reduction by a factor of 200, which brought the life time search 

down to a little less than a year, but this amount of time was still not 

feasible for a computer search. Golay remarked in a footnote that a 

complementary pair might be thought of as being composed of two inter- 

9 

laced half length codes, in this case n=13. Since it was obvious the 
total number of ones in both of the interlaced pair must equal the number 
of ones in one of the complementary code pair, equation 2. 5 was 
applied again, this time for n=13. The unordered weight possibilities 
for n=l3 were (9, 7), (9, 6), (6,4), (7,4). In order to be compatible with 
the (16, 15) pair for n=26, the (9, 7) pair was selected for the code con- 
taining 16 ones and the (9, 6) pair was chosen for the other code which 
contained 15 ones. 

The basic problem had now been reduced to all possible combinations 
of the four thirteen bit sequence families. This was equal to 

^ ^ x ■ - ~ 715x715x1716x1716 or approximately 

1.5x10^^ possibilities. Although this reduced the problem by one more 
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order of magnitude it was still much too large for a practical search. 

Dr. Golay in a private communication with the author pointed out the 

possibility of using the supplementary characteristic as an additional 
14 

screen. Theorem 6. 1 gives the necessary and sufficient condition 
that complementary codes sequence quadruple form be supplementary 
sequence quadruples. The codes of length 13 were therefore categor- 
ized according to their first 3 and last 3 bits. The number of like 
pairs at spacing twelve, at spacing eleven, and at spacing ten was com- 
puted for each of the 13 bit sequences. For example, 0111100111011 
has 0 likes for spacing twelve, 1 like for spacing eleven, and 3 likes 
for spacing ten. The block number given to 0111100111011 was there- 
fore 013. Out of sixty four possible combinations only eighteen different 
block numbers were generated for each of the twrf number sequences of 
all possibilities of thirteen bits with nine ones and with seven ones. An 
examination of Figure 7. 3 reveals that there are a widely varying number 
of members within each of the blocks. The supplementary characteristic 
states that the total number of like pairs for each spacing must equal the 
total number of unlike pairs for the same spacing. With four codes at a 
spacing of twelve there are 4 possibilities, therefore two of them must 
be likes; at a spacing of eleven there are eight possibilities, therefore 
four of them must be likes; and likewise at a spacing of ten there are 
twelve possibilities of which six must be likes to satisfy the supplemen- 
tary property. Block numbers were added for each possible combination 
to see if the total were 246, and in cases where the total was 246 those 
particular codes within these blocks were then sent to the next necessary 
but not sufficient condition for the codes to be complementary. If the 
total was not 246, the blocks of codes were rejected. The power of this 
screen is shown in the example where say two of the blocks had code 
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I 



I 






^ f 




I 




numbers 000, 000. The only possible combination to add to 246 

would be 123, 123 for the other two codes. This cuts the total number 

of possibilities in this particular block combination down from 
10 ^ 

1716x1716 = 2. 9x to only 64 possibilities. Although this admittedly 
was the most extreme example, this screen reduced the number of 

9 

possible codes checked to about 2x10 or further reduction of three 
orders of magnitude. 

The next necessary condition to be checked was the parity test, 
equation 2.4. All code quadruples which satisfied the parity test were 
then sent to the necessary and sufficient like pair, unlike pair check for 
each spacing. The like, unlike check subroutine was quite long and also 
involved word unpacking. This subroutine which was ignored in the 
calculation of run time increased the computer search time from a calcu- 
lated 70 hours to an actual run time of the order of 75 hours. 

The half length codes in the program were interlaced in the same 
manner as was done in the operations group formulation. This method 
of combination made the parity check easy to calculate since the bit 
positions in the computer words were the same as the sequence posi- 
tions of the codes. Modulo 2 addition was directly applied to the code 
quadruple to form a vector of all ones, providing the parity test held. 

If the code vector was not all ones the code quadruple was rejected. 

It is to be noted that by selecting the ordered (16,15) pair for the 
number of ones in the A and B code, the operations Ci, Cs , and C 
were eliminated from possible redundant consideration. These were 
eliminated by their requirement for (10,15), (16,11), and (10,11) 
ordered pairs of ones respectively, based on equation 2.5. The inter- 
lacing of the half length codes with ordered pairs of (9, 7) ones and 
(9,6) ones respectively eliminates from consideration the operations 
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Flow Chart for Computer Search 
For Kernels of Length 26, 

FIGURE 7. 1 
93 



Ti , Ts , T and Ai since these would ha,ve required (7,9») (9,6), 
(9,V)(6,9), (7,9)(6,9) and (4,7)(4,6) ones respectively. The only basic 
transformation not deleted by controlling the weight in each of the code 
quadruples was As. This is seen to be true because if the pair 
(9,7){9,6) has the operation As performed on it, the result is (9,6)(9,7). 
If this code is now transformed by the trivial exchange operation, E, a 
code pair of the original form is generated. The scheme which was used 
to eliminate the As redundancy will be discussed late in this chapter. 

Figure 7. 1 is a rough flow diagram for the program used in the 
exhaustive search for kernels of length 26. With modifications in data 
as necessary, this program was also used for exhaustive searches of 
composite codes of 16 and ^20. 

A brief description of the blocks in Figure 7. 1 will give the reader 
some familiarity with the programming philosophy used in this search. 

The actual programs as written for the CDC 1604 are given in Appendix V. 

The number generation subroutine is able to generate any number up 
to 48 bits in length (word size of the 1604) with any number of ones up 
to a maximum of 48. This is accomplished in the 1604 by inserting a 
number of the form 00000. . . 00001 1 1 1. 1 1 1 into the accumulator, 

where the total length of the number is 4 8 bits and it contains p ones 
and has a code length of n. This number is then shifted left 48-n bits. 
This shifted result is then checked to see if it is negative, which indicates 
a one in the most significant bit position. A procedure is now set up to 
count a one or zero and then shift left one bit and repeat the process for 
the entire n bits. If p ones are counted the number is stored, if not it 
is rejected. A number one larger than its predecessor is inserted into 
the system and the process is repeated until the ones are all in the 
uppermost bit positions. 
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Blocks Z,3,4 are used to compute a supplementary characteristic 
number for each of the 715 thirteen bit numbers with 9 ones and also 
for each of the 1716 thirteen bit code numbers with 7 ones; then in 
accordance with their supplementary numbers these are stored into 
group blocks. The supplementary number was derived by using a mask 
to expose for consideration only the first 3 and last 3 bits of each of 
the 13 bit code numbers. These masked numbers were then compared 
with an exhaustive list of all possible combinations of zeros and ones in 
these 6 bits. These 64 possible combinations were in an M to one 
correspondence with a list of 18 possible supplementary characteristic 
numbers (see Figure 7.2). This correspondence allowed the temporary 
attachment of this supplementary number to each of the codes. The code 
numbers were rearranged into blocks according to supplementary number 
and at the same time had their temporarily attached supplementary 
number deleted. This part of the program is purely for putting the data 
into a useable state and is used only once per run. The rest of the pro- 
gram is highly iterative. 

Consideration of the (I II III IV) form of the possible 26 bit 
codes, shows that each of the 18 blocks of I, II, III, and IV must be 
compared against each other indicating 18x18x18x18^105,000 possible 
block comparisons. This was cut down somewhat by use of a slightly 
different technique. This was accomplished by loading the computer 
accumulator with 246 and from this subtracting the first block number 
of I, then subtracting the first block number of III, and then subtrac- 
ting the first block number of II. Rather than subtracting the Block 
number of IV, the difference obtained was instead checked against all 
all possible block numbers of IV to see if it was listed. Since if the 

block number did exist only one could exist. This procedure cut the 
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possible number of block comparisons by a factor of eighteen and left 
only 5800 combinations to examine rather than 105,000, If the block 
did exist all the codes within this (I II III IV) block grouping were 
then given the parity test. If the IV block number did not exist II was 
stepped ahead one block number and the process repeated. After II 
had cycled through its 18 block numbers, III was then stepped ahead 
one block number and the entire process repeated. Similarly when III 
had cycled its 18 block numbers, I was stepped ahead and the process 
repeated. To avoid the As operation redundancy, whenever I was 
stepped up one block number III was started from this same block 
number rather than at the first block number. This avoided the block 
number combinations (I III) ^ and (III I) both appearing except when 
block I was equal to block III. This block number cycling allowed a 
vantage breakpoint for partitioning the program into suitable size in- 
crements for computer run times. 

The parity check was made by adding(Modulo 2) a code from I 
to a code from III to a code from II. This resultant was then checked 
to see if it existed as a code in the list from the IV block. If it did 
exist it was complemented and then sent to be unpacked for the necessary 
and sufficient like, unlike test. If it did not exist another set of codes 
was sent in and the process repeated until all possible code combinations 
in (I II III) were exhausted. 

The methods used in the unpacking and like- unlike subroutines are 
quite straight forward and will not be amplified here, although it is to be 
noted that the like- unlike subroutine is not limited in the lengths of the 
codes it can test, up to the machine storage size, whereas all the other 
subroutines, since they use packed words, do have quite restrictive 
code length limitations. 
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This program was first used to make an exhaustive search for 
codes of length 10 for test purposes since all kernels of this length 
were known. After this satisfactory checkout the exhaustive search for 
codes of length 2 6 was undertaken. 

During the first part of the computer run for n=26 Dr. Golay in- 
formed the writer by personal correspondence that he had discovered a 
kernel of length 26, using a ”by-hand” technique. The exhaustive 
search revealed only a transformation of the code pair that Dr. Golay 
had discovered. Therefore it was proved that only one complementary 
kernel of length 26 existed and that it was {ignoring allowable trans- 
formations) 

A=0100 1101 11101011 110011 1010 
B=101 100100001 1111 1100111010 . 
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Supplementary 

Characteristic 

Number 

0 

20 
1 10 
1 

1 1 
21 
101 
111 
121 
2 

12 

22 

102 

112 

122 

13 

103 

123 

Total 



Number of 13 
bit codes with 
nine ones 


Number of 13 
bit codes with 
seven ones 


14 


70 


14 


70 


28 


140 


44 


1 12 


44 


1 12 


44 


1 12 


22 


56 


88 


224 


22 


56 


14 


70 


98 


196 


14 


70 


28 


140 


70 


56 


70 


56 


44 


1 12 


22 


56 


35 


8 


715 


1 716 



Number of codes within each 
supplementary block for n=26 



Figure 7. 2 
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CHAPTER VIII 



A PARTIAL SEARCH FOR KERNELS OF LENGTH 34 
The next possible kernels for investigation after length Z6 were 
those of length 34. As will be shown in this chapter an exhaustive 
search for kernels of length 34 was not feasible. However, an impor- 
tant subset of possible kernels was exhaustively searched and a scheme 
which can be adapted to the general search was programmed. No ker- 
nels were found in this partial search. To use the same method of 
attack as that used on n=Z6 did not seem too feasible, both because of 
the length of time involved in the search and because there would be no 
guide to searching, beyond random picking of possibilities until all 
possibilities were exhausted. The length of time involved for an ex- 
haustive search of n=34 was^ estimated fairly reasonably by checking 
the number of possible codes of length Z6 against the number of pos- 
sible codes of length 34, since the time for the Z6 search was known 

with reasonable accuracy. For n=Z6 there were 71 5x715x171 6x1716^^ 
IZ 

1.5x10 code possibilities at one point in the screening process. At 

the same point in the screening process for n=34 codes there were 

1 Z3 76xlZ3 76x1 9448x1 9448 5o 8x1 0^ ^ possibilities, which is a ratio of 
4 

4x10 to one as compared to n=Z6, A time estimate from the Z6 

4 6 

length code to the 34 length code was therefore 75x4x10 =3x10 hours. 
This length of time was of course not practical for a search of an ex- 
haustive nature. 

The use of the cyclic complementary property as a screen was 

pointed out by Dr. Golay as highly practical because some reasonable 

19 

guess could be made as to which cyclic sets might contain a kernel. 

This would lead to a hybrid ”by-hand” and computer search. It would 

also be easier to document the areas searched by the cyclic method as 

compared to the supplementary method and would therefore avoid the 
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duplication of effort by future n=34 length searcherSe As will be dis- 
cussed at some length, there is in addition an even more important 
reason leading to the decision to substitute the cyclic complementary 
property for the supplementary property as a screen. 

The generation of cyclic sets which are to be utilized in the search 
for new kernels, is accomplished in a somewhat unusual fashion. This 
method of classifying, puts within the same set all possible numbers 
which are cyclic permutations of one another, as is to be expected; 
however, there are also contained within the same set all numbers 
which are formed by the removal of every other bit in a cyclic fashion 
until all bits are used. As an example the five bit number 11010 will 
be used to show all the numbers contained within its set. First, 
selecting every other bit, starting with the first bit, and repeating on 

the result gives: 1. 11010 

2 . 10011 

3. 10101 

4. lllOO 

1 . 11010 

and second, by permuting the bits of each of these in the normal cyclic 
permutation manner yields 

1. 11010, OllOl, 10110, 01011, 10101 

2. 10011 , 11001, 1 1100, OHIO, 00111. 

It is not necessary to generate the cyclic permutations from 3 and 4 

since they are contained within 1 and 2. 

Counting cyclically the number of likes for each spacing gives the 
code number 1331 for 1 and the code number 3113 for 2. These 
code numbers start with a spacing of 1 and include up through a spacing 
of n- 1 or four in this particular examplco The code number 3113 
means that at a spacing of one there were 3 likes in 10011 counted 
cyclically, at a spacing of two there was 1 like, at a spacing of three 
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there was 1 like, and at a spacing of four there were 3 likes. The 
example just given is somewhat trivial since all binary numbers of 
length 5 with 3 ones fall into the same set. However, this is not the 
case for n=26 as there are 7 sets for codes of length 13 with 9 
ones and 76 sets for codes of length 13 with 7 ones. 

An examination of the kernel of length 26, shows that its sequence 
quadruple of codes came not from four different cyclic sets as one 
would expect but came instead from just two cyclic sets. 

Using the kernel form 

A=01001 101 11 1010111100111010 

B = 10 1 1 001 00001 1 1 1 1 1 1 00 1 1 1 0 10 one can decompose this into 

1=0010111110111 
11=0010110011101 
III= 1 1010011101 11 
IV= 0010111 100010 

Using I to form its cyclic set gives: 

1 . 0010111110111 
2 . 0111111001101 

3. 0111011111010 

4. 0101100111111 

1 . 0010111110111 

Now take 3 and write it twice, this gives 01110111(1101001110111) 

11010. Picking off from the ninth to the twentieth second bit gives a 

code which is the same as III above , Therefore I and III are 

members of the same cyclic set. Similarly take II and form its 

cyclic set; this gives: 

1 . 0010110011101 
2 . 0110111001010 

3. 0111000101011 

4. 0100001110111 

5. 0001111100101 
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6 . 0011011011100 

7. 0101110011010 

8 . 0010100111011 
9 . 0110101000111 

1 0 . 0111011100001 
1 1 . 0101001111100 
1 2 . 00011 101 10110 
1 . 00101 10011101 

To make a comparison between II and IV the complement of IV will 
have to be used in order that it will have 7 ones as does II« 

IV= 1101000011101 

Writing 4 from this list of 12 twice gives 0100001 1101(1101000011 101)11. 
Picking off from the twelfth to the twenty fifth bit gives IV, which in- 
dicates that II and IV are hrom the same cyclic set. 

A check of kernels of length 10 to see if their quadruples came 
from cyclic sets in pairs, as the 2 6 kernel did, proved that both ker- 
nels of length 10 did in fact come in pairs. This was meaningless 
however, since as was seen in the example of length five with three 
ones, there was only one cyclic set possible for each pair. 

In the hope that if the kernel existed in the n=34 case, it would 
also be formed with I, III from one cyclic set and II, IV from another 
cyclic set, the program was modified. This change along with giving an 
exhaustive search for code possibilities taken two at a time from cyclic 
sets would also give an estimate of the computer run time for the ex- 
haustive search in general using the cyclic complementary property. 

This was the primary reason for the change in the program. 

Eight octal digits were used to designate the cyclic complementary 
number for each of the 12,376 possible half codes of length 17 with 
11 ones and for each of the 19j.448 possible half codes of length 17 
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Flow Chart for 

Cyclic Complementary Designation. 
FIGURE 8. 1 
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with 10 ones. Only eight digits were necessary as each cyclic count 
of spacing i actually counts both the i spacing and the n+l-i spacing 
simultaneously. This was demonstrated by step 8 in the proof of 
Theorem 6.2. The example of the code of length five given earlier in 
this chapter showed this symmetric property, as the code numbers were 
1331 and 3113. The first two digits of these cyclic complementary 
designators contain all the information available. 

The revision made to the computer program was not too extensive 
as only blocks 2, 3,4, 5 in Figure 7.1 were modified to handle the 
cyclic complementary property rather than the supplementary property. 
Of course all blocks were modified to handle codes of length 34 rather 
than codes of length 26. 

One difficulty encountered was the lack of computer memory for the 
large blocks of data, and it was necessary to store the cyclic sets on 
magnetic tape rather than in the main memory core as the 26 case was 
handled. The cyclic sets were then called in two at a time for checking. 
This slowed the search down a little, but by no more than 1% of the 
total time. Figure 8. 1 is a block diagram of the process of obtaining 
the cyclic complementary numbers and it will be worthwhile to remember 
the memory size restriction when reviewing the procedure used. 

A general description of Figure 8. 1 will show the differences in 
the method used in the search for n=34 kernels as compared to that 
used for n=26 kernels. The procedure in the n=34 data formation was 
identical with that of n=26 as far as the initial possible number genera- 
tion was concerned. (Block 1 in Figures 7. 1 and 8. 1). There was 
an immediate departure from the old method in blocks 2,3, and 4 
where one of the 17 bit numbers with 11 ones was shifted left 17 bits 
and then added to itself. This formed a 34 bit number which was 
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actually the 17 bit number written twice in time sequence. A mask 
was then put over the last 17 bits and this masked number was used, 
with an equality search of the 12,3 76 generated 17 bit numbers, to 
identify this particular number in the list. After the number equal 
to the reference number was located in the generated list it was tagged 
a one, indicating cyclic set one. The whole 34 bit number was then 
shifted right one bit and a mask used to again pick off the last 17 bits. 
The search through the list and tagging with the group number was re- 
peated. After cycling through all 17 possible cyclic codes and tagging 
each, the original 34 bit number was now repeatedly long right shifted 
one bit and then right shifted one bit, A long right shift saves the bit 
which is pushed off the A register and stores it in the Q register, 
while the right shift just pushes the extra bit off the end of the A reg- 
ister, Therefore every other bit of the 34 code was saved in Q and 
this formed a new 17 bit code within the same complementary cyclic 
set. This code was then formed into a double length code and the list 
of all possible 17 bit numbers with 11 ones was searched and tagged 
with a one and the process which was first described to generate cyclic 
permutations was repeated. This process of first forming 17 cyclic 
permutations of a code and then taking every other bit, was continued 
until the possible code which was originally operated upon reappeared. 
Then cyclic set 2 was started and tagged accordingly, and so on 
through 98 sets. The 19 1 448 codes of length 17 with 7 ones were 
classified in exactly the same manner and formed 150 cyclic comple- 
mentary possibility sets. 

After each series of sets were formed a mask was used to search 

the set tag numbers and pick off each set for transfer to magnetic tape. 

To insure that the pick-off was exhaustive, a zero was inserted on the 
list in place of each removed tagged number, 
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The sets were now ready to be called back into the computer for use 
in the search. Normally only one set was called in at a time because 
the pair from the I, III set was checked against all the II, IV sets 
before a new I, III set was called into the computer. The I, III sets 
contained the 11 weight codes while the II, IV set contained the 7 
weight codes. When a set was called into the computer each cyclic 
permutation subset of 17 possible codes was given a cyclic comple- 
mentary number to categorize the entire subset. Each group contained 
either 4 or 8 of these cyclic subsets. These subset numbers were 
formed by counting the likes for the I, III codes and the unlikes for the 
II, IV codes for the eight unique spacings. The codes within the subsets 
were checked against each other only if the likes equaled the unlikes of 
the cyclic complementary subset numbers for all spacing, or Lj + 



hii = Un + 



u 



IV 



Figure 8, 1 is joined to Figure 7, 1 just beyond block 5 and takes 
the place of blocks 1,2, 3, 4, 5 in Figure 7. 1, The same feedback paths 
shown in Figure 7. I to insure the checking of all codes within a set 
and to feed in a new set after a check has been completed, are still in 
operation for the same purpose. Appendix VI contains the computer 
program for Figure 8, 1, 

The total run time for the n=34 search, taking I, III from the 

same cyclic set and II, IV from another cyclic set was approximately 

15 hours. No codes were found. An extrapolation from this run time 

5 

to the exhaustive run time for n = 34 gives 1 5x 1 50x98 2, 2x 10 hours, 

which although one order of magnitude less than the supplementary 
method, is still not reasonable, 

A second partial search for n=34 kernels was attempted after 

observing that the kernel of length 1 0 
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A=0101000011 
B=00001001 10 

with its center quad removed yields 

A^=01010011 
B^=000001 10 

which is a complementary pair of length 8„ Breaking A^, B^ into its 

(I II III IV) configuration gives 

1=0001 

11=1011 

111=0001 

IV=0100 

where (I, II) and (Ills IV) are both complementary pairs. 

An attempt was made to draw an exact parallel from this observa- 

tion to codes of length 34. One possible solution to equation 2.5 which 

determines the number of ones necessary in each of the kernels of 

length 34 is (21,18). If this same equation is applied to codes of 

length 17, two of the possible solutions are (11,10) and (11,7). The 

(11, 10) pair when interlaced will total 21 ones, and similarly the 

(11, 7) pair when interlaced will total 18 ones. Whereas the kernel of 

length 10 with the quad removed had a sequence quadruple form with 

codes of length 4, the 34 code with the quad removed would have a 

sequence quadruple form with codes of length 16. Appendix II gives 

a list of all possible codes of length 16 with weight 10 in the A code 

and weight 6 in the B code. This list does not include the results of 

the time reversal transformations which have to be included since weight 

is invariant under T^ , Tg and T. It was observed that if these codes 

of length 16 were interlaced to form codes of length 32 and then sep- 

erated in the middle and the quad A=ll, B = 10 inserted, all conditions 

as far as the number of ones required in both interlace codes of 17 to 
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form codes of length 34 would be satisfied. An example might serve 

to clarify the preceding statements. 

0=0101111111000110 

D=1001001100001010 

is a complementary pair of length 16^ Using the interlace method to 

form a pair of length 3 2 gives 

A=01 100011 101011 11 1010000001101100 
B=001 101 10111 1 10101 11 10101001 1 1001 . 

Splitting this A,B pair in the center and adding the quad to this center 

give s 

A^=01 10001110101111 11 1010000001101100 

B^=001 101 101 1111010 10 1111010100111001 . 

Breaking this pair into standard (I II III IV) form gives 

1=01011111111000110 

11=01010000111001001 

111=0101 11 111 11000110 
IV=10101111000110110 . 

I has 11 ones. Interlacing I with II, which has 7 ones, yields 
a total of 18 ones. Similarly III has 11 ones, and IV has 10 ones. 
When these are interlaced the code of length 34 has 21 ones. In the 
actual computer run the codes of length 16 with 6 ones were converted, 
by complementing, into codes with 10 oneSo Therefore both the A and 
B codes contained 10 ones« An exhaustive list of these and their time 
reversals was generated. Each code of length 16 was then split in the 
center and a zero inserted for the ninth bit, forming an over all code of 
length 17 with 10 ones. A second list was then formed, using a one 
rather than a zero for the inserted ninth bit, forming codes of length 17 
with 11 ones. 

These two lists of numbers were then fed into the program of 
Figure 7. 1 at block two, and the regular supplementary check program 



was used. The codes with 11 ones were used for (I, III) and the 
codes with 10 ones were used for (II, IV)o When these codes were 
combined in the interlace scheme, I interlaced with II gave 21 ones 
and III interlaced with IV gave 18 ones. All four interlace combina- 
tions (I II m IV), (II I III TV), (I II Tv III) and (II I Tv III) were 
run exhaustively. The run time was approximately three hours and no 
codes were found. These partial searches for n=34 are documented 
here because they represent two fairly obvious approaches to the prob- 
lem, and it would be very wasteful of time for someone to duplicate this 
effort. The program as given in the appendix can with a slight modifica- 
tion be used for an exhaustive search for kernels of length 34 or for a 
partial search based upon other ^judicious guesses**. 
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CHAPTER DC 
CONCLUSIONS 



The study which has been presented here had two principal objec- 
tives. The first objective was to make an exhaustive search for kernels 
of length 26 to determine whether they existed, and if they did exist, 
to determine the number of such kernels. The second objective was to 
develope a methodical scheme of code decomposition from composite 
codes back to their generating kernels. We shall summarize in this 
chapter the ways in which each of the major portions of this dissertation 
is related to the accomplishment of these objectives. Also a number of 
suggestions for further research will be presented. 

The operations group and its family of theorems were originated to 
formalize the elimination of possible redundancies in the search for new 
kernels. The predetermination of the number of ones in each of the 
sequence quadruple vectors was a vital screen that eliminated many 
such redundancies. Also, the starting of both the I and III vectors 
from the same possible supplementary code groups eliminated an ad- 
ditional half of the possible codes. The theorems in the operations 
group chapter also were useful in the code decomposition part of the 
problem, since they show the necessary and sufficient conditions for 
the (I, II) and (III, IV) pairs to be complementary. Theorem 3.4, 
which states that if the sequence quadruple form is (I II I nj then 
(I, II) must be a complementary pair, is most helpful since the process 
of checking for complementarity is very tedious by hand methods. 

Hamming distances were the next major topic considered. They 
were found to be useful in a number of the proofs in the Hamming vector 
chapter. They were also a great aid in checking hand decomposition, 
since it was easy to determine whether D(I, II)=D(III, IV)=n/4 for compo- 
site codes or whether they had the characteristic distance in the case of 
kernels. Ill 



Hamming vectors were the most powerful tool developed in this 
paper for code decomposition^ The theorems cover all known possible 
decompositions and also give strong support to the conjecture that ker- 
nels of length 50 do not exists The Hamming vector concept should 
allow a fuller understanding of complementary sequences in general, 
due to its characteristic qualities such as anti- symmetry, and also be- 
cause there are only four possible configurations of Hamming vectors 
for any code pair. Similarly in sequence quadruple form it is easy to 
check whether a decomposed code may be complementary by using its 
Hamming vector. 

Several exhaustive searches for codes of various lengths were con- 
ducted during this research. The most important of these was the 
search for kernels of length 26. This search revealed that only one 
kernel existed. An exhaustive search for codes of length 16 and 20 
revealed that all codes of these lengths were generated from shorter 
codes by standard methods. This gives considerable strength to the 
conjecture that these generating methods are the only ones. 

The possibilities for future research in this field of complementary 
sequences are many. Some of these follow directly from this paper 
while others lie quite far afield from the ideas exposed here. An inter- 
esting extension of the present work would be an exhaustive search for 
kernels of length 34. It is felt that the cyclic complementary screen 
which was used in the partial search is not powerful enough, and a 
better screen is needed before this search can be conducted on an ex- 
haustive basis using today’^s computers. The writer had considered the 
combination of both the cyclic and supplementary properties as screens 
in the same search. However, a full study was not made of the gain 



which might be obtained by this screening procedure. The programming 
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of such a screen would not prove too difficult and, since the cyclic 
property would logically be the first screen of the two to be used, the 
areas of search could still be catalogued quite easily. 

If it would be possible to determine, before the search, the charac- 
teristic Hamming distances D(I, II) and D{U^ IV) of kernels, these 
would act as a most potent screen in the search for new kernels. It is 
therefore felt that further research on characteristic Hamming dis- 
tances of kernels could prove quite fruitful for future applications. 

It was conjectured that kernels of length 50, 130 and 338 do not 
exist, based on several pieces of evidence, which although strong are 
not conclusive. Further investigation along these lines might offer a 
conclusive proof without the necessity for an exhaustive search, since 
an exhaustive search for these length codes approaches the impossible. 
During the investigation of kernels of length 50 the writer noticed an 
oddity which might be worthy of further investigation. 

One kernel of length 10 is 
A=1001010001 
B=10000001 10. 

Divide both the A and B codes into segments of length 2, then 
A=10 01 01 00 01 

B=10 00 00 01 10. 

Now let the symbols 10=A, 01=A, 11 = B, and 00=B. Forming a code 

of length 50 by this method of using symbols of portions of the quad 
to represent kernels of length 10 gives 

Ki=aXaba 

K2=ABBAA. 

This code pair of length 50 satisfies the parity check, and also equa- 
tion 2. 5 for the number of ones in both full codes and in the four half 
codes. Several examples were checked in this fashion for both kernels 
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of length 10, and all satisfied the same conditions as did this example, 
but of course none of them satisfied the necessary condition for count- 
plementarity. Along with helping to search for codes of length 50, this 
offers a possibility for research in that all kernels seem to have 
various characteristics in common^ However, the quad seems more 
flexible. Of the four known kernels, the quad is the only one thus far 
discovered that can be used to build up composite codes by one method 
of generation and then these codes can be decomposed by a different 
method. Therefore a study of kernel characteristics might lead to 
easier ways of generating new kernels, or at least of determining 
whether they exist for the various possible lengths. 

Another possible field of research is cyclic complementary codes. 

■s 

A communications system is more likely to use a cyclic complementary 
carrier modulation instead of straight complementary modulation due to 
the continuous nature of the carrier. Therefore an investigation of 
these codes seems quite in order. 

Two other topics worthy of research in the complementary sequence 
field are the various correlation functions of complementary sequences 
and the frequency spectra of these sequences. For secure communica- 
tions purposes it is important to have a code that is noise-like. A 
uniform distribution of ones and zeros gives a noise- like appearance to 
the code. For example, the complementary pair 

A=11010001 

B=11011110 

is a typical complementary code pair of length 8. If the bits of the A 
code are considered in a cyclic fashion one, two, and three bits at a 
time, the results are remarkably uniform as is seen in Table 9. 1. 

This uniformity in output would make it extremely difficult for someone 



monitoring the code to decide if he had noise or an actual signal on his 
receiver. However, it is seen that the B code is not nearly so uniform 
and would be easier to detect as a signal. 



A code 



Bits 


Number 


Two Bits 


Number 


Three Bits 


Numbe 


1 


4 


1 1 


2 


1 1 1 


1 


0 


4 


10 


2 


1 10 


1 






01 


2 


101 


1 






00 


2 


Oil 


1 










100 


1 










010 


1 










001 


1 










000 


1 


B code 




* 








1 


6 


1 1 


4 


111 


2 


0 


2 


10 


2 


1 10 


2 






01 


2 


101 


2 






00 


0 


Oil 


2 










100 


0 










010 


0 










001 


0 










000 


0 



Noise uniformity of a complementary 
pair of length 8, 

TABLE 9e 1 

An investigation of a large number of codes by the above technique would 
be useful in determining the most nois e-like code pairs, where noise- 
like is defined as a uniform output with bits taken one, two, three, etc. 
at a time. This might in turn lead to some general characteristic of 
most-noise- like complementary sequence pairs, 
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Another important study would be the autocorrelation function of 
each code pair. This is synonymous with the output of a code detected 
by its matched filter. For example, the autocorrelation of the code A 
given above is 

1 0 -3 0 -1 0 - 1 8 -1 0 -1 0 -3 0 1 

where the outputs in time are read from left to right. The autocorrela- 
tion of the B code is 

- 1 0 3 0 1 0 1 8 1 0 1 0 3 0 -1. 

The autocorrelation functions are important from two standpoints. First, 
if a portion of the system should need to be shut down for repairs or 
routine maintenance, a good autocorrelation function would allow opera- 
tion on just one carrier. Second, the frequency spectrum of the trans- 
mitted signal is the Fourier transform of this autocorrelation function. 
Therefore, one way to study the spectrum is first to autocorrelate and 
then to take the Fourier transform of the transmitted signals. The more 
widespread the spectrum, the more difficult is a jamming procedure, 

A study of the crosscorrelation functions of complementary sequence 
pairs is also very important because of the difficulty of isolating the A 
code RF carrier from the B code receiver. It was pointed out in 
chapter 4 that because the A code and the B code were orthogonal, at 
exact match cross talk would be no problem. However, when the signals 
are not at exact match cross talk might be a large source of noise. An 
inve stigation of the c ros scorrelation functions of the A and B code 
pairs might lead to a class of complementary sequences where the 
problem of crosstalk is minimized. 

Another class of correlation functions which are important are 
those where cyclic errors are created in the received signal detected by 
the matched filter. These errors could be caused by a linear shift in 
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the phase of the received signah There may be certain classes of 
complementary sequences v/here the loss in coding signal-to-noise ratio 
is minimized when noise is caused by an error of this type* This par- 
ticular class of signals would be extremely important in an air search 
radar application* 

A study of the maximum number of different composite codes for 
any length would be an important contribution to the field, because for 
anti- jam reasons it is very important to be able to change codes, and 
the more codes with good characteristics that are available the better. 

The last study to be recommended is an analysis of the autocorrela- 
tion functions, frequency spectra, cross correlation functions , and 
linear phase shift correlatipn functions for codes of lengths from 100 
to 200. The purpose of this investigation would be to determine if 
certain kernels are capable of forming composite codes of better char- 
acteristics than other kernels for radar applications. 

In conclusion, the writer feels that complementary sequence pairs 
will prove very important in future applications to both radar and com- 
munications schemes^ However, much research must be done to pick 
optimal classes of these codes for such applications. This paper has 
offered many tools to be used by future investigators of complementary 
sequences. These tools have found important applications in their 
present state, but additional investigations along the lines of operations 
groups, Hamming weights and Hamming vectors should further in- 
crease their usefulness. 
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GLOSSARY OF SYMBOLS AND TERMS 



(A.B) 


A complementary pair of binary sequences, 
see equation 2^ 2. 


(I II III IV) 


A complementary pair of binary sequences in 
sequence quadruple form, see equation 3« 1. 


D(U,V) 


The Hamming distance of two binary vectors , 
see equation 4« 1„ 


H(U.V) 


The Hamming vector of two binary vectors, 
see the first page of Chapter 5« 




Means interlace two Hamming vectors in a pre- 
scribed manne r, see Theorem 5.11, step 2. 

Is transformed into. 




The time inverse of a binary sequence or 
vector A. 


Alte r 


Means to complement every other bit of both 
codes of a complementary code pair. 


Complementary 

code 


A pair of binary sequences which satisfy 
equatipn 2. 2, also called a complementary 
sequence pair, a complementary sequence 
and sometime just a pair or code. 


Complement a 
code 


Means change the sign or take the complement 
of a binary sequence. 


Composite code 


A complementary code which is reducible to a 
shorter code Length by standard methods. 


Cyclic comple- 
mentary code 


A pair of binary sequence which satisfy 
equation 6. 3. 


Hamming Distance 


The number of bits in which two binary sequence 
differ, notation is D(A,B), see equation 4. 1. 


Hamming vector 


A binary vector or sequence formed by the 
modulo 2 sum of binary sequences. 


Hamming weight 


Hamming distance of a binary sequence or vector 
from the null vector. It is the number of ones 
in a binary vector. 


Ke rnel 


A set of complementary codes which is irreducible. 
The term is sometimes used for a member of the 
set. 


Quad 


A kernel of length 2. 


Sequence 

Quadruple 

Supple me ntary 
code 


A complementary code pair expressed as in 
equation 3. 1 ® 

A quadruple of binary sequences which satisfy 
equation 6. 1 . 


Time inverse 
Weight 


Means to make the first bit of a binary sequence the 
last bit, the second bit next to last, and so on; also 
called time reverse. 

See Hamming weight. 
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EXPLANATION OF FORMAT OF THE APPENDICIES 



Appendix I on matched filters is self-explanatory. Appendix II 

is the result of an exhaustive search for codes of lengh 16, There were 

96 code pairs found by the computer after the program screens had 

eliminated most of the redundancies. The format of Appendix II gives 

the A code, the B code, H(A, B), and H[Ti (A,B)] which signifies 

time reversing the A code. It is worthwhile to note that each code pair 

and the pair formed by the time reverse of A were decomposable. 

From these 192 possible decompositions 48 were interlace, 48 were 

r 

time sequence, and 96 were 2 special method, with 48 each of two 
different kinds. 

Appendix III is the re^sult of an exhaustive search for codes of 
length 20. The format of Appendix III gives the A code, the B code, 
if it is required, the time reversal operation used to put H( A, B) in 
standard form, H(A,B) and from which of the 2 kernels of length 10 
the generating codes came. All of the 24 code pairs or their time 
inverse were decomposable. There are 12 from each kernel, 8 of 
each of these interlace and the other 4 are time sequence. 

Appendix IV is a diagrm of all the subgroups of the operations 
group. It starts with the subgroup of order 1 and works up to 
those of order 16. An x in the row indicates that the operation at 
the top of the column is a member of the subgroup. 

Appendix V is the CDC 1604 computer program in AR format 
and also in machine language for the 26 length search. There are 
slight modifications in the various sections of the program required 
for other length codes. Appendix VI is also in both AR and machine 
language, and is for codes of length 34. This program parallels 

Figure 8. 1, whereas the program in Appendix V parallels Fig. 7. 1. 
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APPENDIX I 



MATCHED FILTERS 

A Matched Filter is by definition a filter which maximizes the peak 
signal- to- noise ratio. Deterioration of the signal wave form is accepted 
in order to obtain the desired maximum ratio. 

One tool needed for the derivation of the matched filter character- 
istics is the Schwarz inequality which is a special case of the Holder 
22 

inequality. One representation of the Schwarz inequality is 




x(w) y (w) dw 



2 Too 

< J |x(w) 

- 00 




y(w) 




(A 1.1) 



It is to be noted that the equality holds when y(w) is the complex 
conjugate of x(w). Rewriting equation A 1. 1 gives 




2 2 

x(w) I dw J |y(w) dw 



(A 1,2) 



The derivation for the matched filter given here closely parallels 
23 

that of Turin. 

Let f(t) be a signal impressed across a filter whose frequency 
characteristic is H(jw), F(jw), the frequency spectrum of f(t), is 
given by the Fourier transform of f(t)o 



F(jw) = 



L «t) « 




G(jw) is the frequency spectrum of the output after passing through 
the filter 



G(jw) = H(jw) F(jw), 
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or transforming back to the time domain the output signal voltage is 



F(jw) e (A 1.3) 

Since at some time, t = A, g(t) must be a maximum, 

r °° 

max g(t) = g(A) = J_^H(jw) F(jw) e (A 1.4) 



To complete the derivation it will be necessary to obtain the total 
noise power, and the total power in the signal. If the noise is assumed 
to be white noise of N watts per cycle, the output noise power density is 
N = N I H(jw) 1^, 

LI "> 

and the output noise power is therefore 



OO J, 00 

J N df = J N iH(jw) I df. 

-.00 o - 1 w / I 

The total energy contained in the signal is 



00 




f^(t) dt, 



(A 1. 5) 



which is also 



r 00 

E = J |F(jw)rdf, (A 1.6) 

by PlanchereTs Theorem. It is to be noted that E is a constant not 
dependent upon the filter used. 

The output power signal- to- noi se ratio for the maximum output is 
the ref o re 
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5 






C 00 

I i F(jw) H(jw) 



N 



N 



O 00 

L |H(j 



jw) l^df 



Dividing this equation by a constant will not effect the time of maximum 
signal output; therefore dividing by E, equation A 1, 6 , gives 



N I g(A) 1^ 
E N 

o 



00 

I F(jw) H(jw) e df I ^ 

f ! H(jw) l^df [ 1 F(jw) 1^ df 

-00 - 00 



(A 1. 7) 



If the right hand side of equation (A 1, 7) is compared *to the 
Schwarz inequality (A 1.2) with x(w) = F(jw) e y(w) = H(jw), 

the two expressions are the same. But for the expression to be max- 
imum y(w) must equal the complex conjugate of x(w), therefore 
H(jw) = F(-jw) e 

A matched filter is therefore a filter whose frequency character- 
istic is the complex conjugate of the signal spectrum to which is matched. 

The output signal for a given signal applied to a matched filter is 
therefore 

00 

g(t)=L |F(jw)!^e _ 

at t = A the output is 

00 

= Joo I • 



This is the same as equation A 1. 6 which is all the energy con- 
tained in the signal. Therefore at exact match all the energy in the 
pulse is in the output signal. 
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APPENDIX II 



Exhaustive Search for Codes 
of Length 1 6. 
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APPENDIX V 



REMARKS ON COMPUTER PROGRAM 

In Appendix V the addresses from 60000 to 60035 are block I 
of Figure 7. 1, All numbers used in the program are octal. The program 
as written is for codes of length Z6, so most of the program deals with 
half length codes of length 13. Words to be changed for different length 
codes are 60000 (mask of ones for half length code as determined by 
weight), 60001 (length of half code minus one [upper half], and length 
of memory word minus length of code [lower half]), 60004 (length of 
half code), 60034 (last sequence in code list). For exactly the same 
reasons respectively change 60015, 60016, 600Z1 and 6003Z. 

Addresses 17000 to 17045 separate the sequences into blocks by 
supplementary number, Ii\dex register Z must be entered with 100 
before starting. Addresses to be changed here are 17000 (mask of first 
3 and last 3 bits), 17005 (contains the total number of sequences for 
both half length codes), 17007 (contains the number of sequences with 
7 ones), 17015 (contains a mask for half length codes), 170Z6 (contains 
the number of sequences with 9 ones), 1703Z, 17033 (have for their 
instruction address the base word 50000 plus the number of half codes 
with 7 ones), 17034 (contains a mask for code length). 

Addresses 10000 to 10106 are the repetitive part of the program. 
The only memory address needing modification here for different length 
codes is 10106 which contains SVN a mask equal to the one half code 
length. 

Addresses 7Z00 to 7Z5Z contain the unpacking subroutine and the 
like-unlike subroutine. 7Z0 1 , 7Z0 7, 7Z1 5 , and 7ZZ3 all contain the 
half code length. Similarly 7Z04, 7Z07, 7ZZ0, 7ZZ3 and 7Z30 all contain 
the full code length or the length modified by one or two to fit the program 



needs. 
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The 15000 to 15100 series is all possible bit arrangements for 
the first 3 and the last 3 bits of codes of length 13, while 16000 to 
16100 is a list of supplementary block numbers in one to one correspon- 
dence with this list. 

The 10300 and 10400 series respectively contain the number of 
codes in each supplementary block for the half length codes of weight 
9 and weight 7, 

16000 to 160Z0 is the final list of the 18 supplementary block 
numbers and is in one to one correspondence with the 10300 and 10400 
series. 
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